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In two dimensional electron systems with Coulomb or dipolar interactions, a direct tran- 
sition, whether first or second order, from a liquid to a crystalline state is forbidden. As 
a result, between these phases there must be other (microemulsion) phases which can be 
viewed as a meso-scale mixture of the liquid and crystalline phases. We investigate the 
transport properties of these new electronic phases and present arguments that they are 
responsible for the various transport anomalies that have been seen in experiments on the 
strongly correlated 2DEG in high mobility semiconductor devices with low electron densities. 

As low density two dimensional (2D) electronic systems with increasingly high mobility have 
become available, there has accumulated experimental evidence of a set of low temperature phe- 
nomena which cannot be understood on the basis of traditional (Fermi liquid based) metal physics 

(More precisely, by low density we mean a 
large ratio between the potential and the kinetic energy in the liquid measured by the parameter 
Vs ~ 1 / n{aBy^ , where n is the areal density of electrons and is the effective Bohr radius.) A 
generic feature of electronic systems is a low temperature evolution as a function of n from a con- 
ducting "Fermi liquid" (FL) state at large n to an insulating "Wigner crystalline" state (WC) at 
low. Much is understood about these limiting cases, so it is possible to make rather firm statements 
about the failure of these states to account for the experimental features. 

Recently, we 13,0, 13; 13; I3| have shown that in ideal 2D systems there must exist new 
electronic micro-emulsion phases at intermediate densities - phases which consist of various sorts 
of micro-phase separated mixtures of regions of WC and regions of FL. This follows from the fact 
that, in the presence of either Coulomb or dipolar interactions, neither a direct first or second order 
transition from a uniform liquid to an electronic crystal is allowed. We have argued that these 
micro-emulsion phases may be responsible for the anomalous behaviors observed in experiment. 

Electronic micro-emulsion phases are distinct thermodynamic phases of matter; transitions 
between them are associated with nonanalytic dependencies of the free energy on T and n. Since 
different micro-emulsion phases exhibit different patterns of phase separation on the mesoscale, 
they also have new forms of hydrodynamics. For example, when the minority phase consists of 



dilute bubbles of FL embedded in a WC host, the result is a relative of a supersolid phase j23l |: 
the hydrodynamics of such a state combines usual hydrodynamics with elasticity theory. The 
hydrodynamic properties of other microemulsion phases are even more complicated both on the 
classical and quantum levels. On the classical level, some of them are similar (but not identical) 
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to classic liquid crystals. 

It is the purpose of the present paper to report a step toward a theory of the low energy, long 
wave-length dynamical properties of these micro-emulsion phases, and to begin the exploration of 
how they are modified in the presence of weak quenched disorder. We have also undertaken to 
compare those aspects of the theory which we do understand with experiments on the 2D electron 
gas (2DEG) in Si-MOSFETs Q, 11 i, |lB,Ji E J1B> and the newer experiments on the 2D hole 



gas (2DHG) in GaAs heterostructures [Ul, UJ, UJ, UJ, lla, ll^, Hi 

Since this paper is long, efforts have been made to make the various sections self-contained. 
Each section begins with a summary of the important results. In Section ^ we discuss some 
of the features of the observed behaviors of the 2DEG and 2DHG at large which we believe 
are incompatible with a Fermi liquid description. In Section ^ and Appendix we review the 
thermodynamic analysis which implies the existence of electronic micro-emulsion phases. Much 
of this discussion is review, although we have expanded considerably upon previous results. The 
central result in this discussion is Eq. [2 (or its more general form in Eq. which expresses 
the free energy density of a general micro-emulsion phase in terms of a non-local surface tension. 
Section HTn is a slight digression, in which the relation is derived between hydrodynamics of a fluid 
in the absence of disorder and the resistivity in the presence of weak disorder. In Section HVl we 
discuss the finite T hydrodynamics of what is in many ways the simplest microemulsion phase, 
a WC bubble fluid consisting of a dilute concentration of WC bubbles in a majority FL sea. In 
Section \V\ and Appendix ^ we analyze the low T limit where the quantum nature of the bubbles 
begins to play an important role. In Section IVII we briefly discuss the hydrodynamics of some of 
the other microemulsion phases. In Section fVIII we consider the case in which there is a bilayer of 
two 2DEG's which permits measurements of the "drag resistance". In Section IVIIIl we compare 
the results obtained in the present paper with corresponding existing experimental observations 
on the 2DEG and 2DIIG. We also discuss what we think are the major shortcomings of some 
of the alternative, Fermi liquid based explanations of the same data that have been discussed in 
the literature. Finally, in Section IIXI we conclude with a discussion of a number of important 
remaining open issues. 

Throughout this article, we work in units in which Boltzman's constant, ks = 1, and the 
Zeeman splitting, gfiB^ = where g and fiB are, respectively, the gyromagnetic ratio and the 
effective Bohr magneton of an electron. 



I. EXPERIMENTAL SIGNATURES OF NON-FERMI LIQUID BEHAVIOR IN 

STRONGLY CORRELATED 2DEGS 

In this section, we summarize some of the experimental results on large 2DEG's and 2DHG's 
that we think are incompatible with Fermi liquid based theory. We mainly focus on the behavior 
of systems with large but also small values of the resistivity p (compared to the quantum 
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of resistivity /i/e^), since in this limit, the predictions of Fermi hquid theory are crisp, and so 

discrepancies with the theory can be documented. (A more specific discussion of some proposed 

FL based theories can be found in Section fVIII CI ) 

To orient the reader, the devices we have in mind have ~ 5 — 20 in Si MOSFETs and 
~ 10 — 40 in p-GaAs. In comparison, the best current theoretical estimates of the critical 

value of Tg at which the energies of a uniform FL and WC are equal is r* ~ 38. 



A. Experiments on single layers. 



a. The metal-insulator transition: Both the 2DEG and the 2DHG exhibit an "apparent zero- 
temperature metal-insulator transition" as a function of the electron density n, or the gate voltage 
in MOSFET's. (See for example. Fig. 1 in Ref. [l| ). Applying a magnetic field, Hp parallel 
to the 2D plane drives the system toward the insulating phase. Thus the critical concentration 
nc(//||) at which the transition occurs increases with Hp The existence of such a transition is, 
perhaps, the most fundamental signature of the breakdown of the Fermi liquid based theory of 
localization, since in this framework there is no metallic phase in the 2DEG at T ^ 0. While it is 
difficult to establish unambiguously that the apparent transition is not actually a rapid crossover, 
the resistivity certainly changes by several orders of magnitude as n varies over a modest range near 
a critical concentration. The range of resistances and temperatures covered by the experiments is, 
in some cases, impressively large. 

Where the resistance p <C h/e'^, were FL theory applicable, perturbation theory in powers of 
1/kpi would be well controlled. Thus it is relatively straightforward to identify behaviors that are 
not consistent with FL theory in the "metallic" regime. (Here, and henceforth, we use the term 
"FL theory" in a somewhat expanded sense to refer to results based on FL theory plus low order 
perturbative corrections of various sorts in powers of l/kpi.) 

b. Temperature dependence of the resistitivity p{T) in strongly correlated (rg ^ Ij metallic ( 
p <^ h/e^) samples: At low temperatures, p{T) is observed to increase with increasing temperature 
by as much as a factor of six in Si-MOSFET's (See Fig. 1 in Ref. P| ) and by a factor of three or 
more in GaAs heterostructures (See Fig. 1 of Ref.|l2j and Fig. 3 of Ref. [l^l)- At T 

the resistance has a maximum at which /i/e^ Pmax ^ piT = 0); then at T > Tmax it decreases 
with T. The resistivity maximum is not terribly pronounced in the data on Si MOSFET's, but 
some of the data taken in p-G aAs exhibit this maximum very clearly, as can be seen in Fig. 1 of 
Ref. 1^ and Fig. 3 of Ref. 13 1. In the latter case p <C /i/e^ in the entire range of the temperatures 
of interest so interference corrections to the conductance are unimportant. 

Such a strong temperature dependence of the resistance is unexpected in the context of Boltzman 
transport in a FL at low T where electron-phonon scattering is negligible. (In many cases, a 
crossover temperature, Tph below which electron-phonon scattering is unimportant can be readily 
identified; for example, in Fig. 3 of Ref. Q|, Tph > Ep is roughly the point at which p{T) 
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has a minimum.) Since the Fermi momentum is much smaller than the reciprocal lattice vector, 
electron-electron scattering conserves the total quasi-momentum and, therefore does not contribute 
directly to the resistance. The resistance is therefore determined by the rate of electron-impurity 
scattering, which typically is a weakly decreasing function of the temperature in a FL. 

c. Parallel field magneto-resistance of strongly correlated (rg ^ 1) metallic ( p h/e^) samples: 
At T < Ef and at sufficiently large r^, both in Si MOSFET's (See Fig. 3 in Ref. jdl and Figs. 1 
in Ref. 0, 0|) and GaAs hetero-j unctions (See Fig. 1 in Ref. ^^), the 2DEG exhibits a strongly 
positive magneto-resistance as a function of i^y, which saturates at H^^ > Ep. 

In sufficiently thin samples, i^y has little effect on the orbital motion of the electrons, so it can 
be viewed as coupling only to the electron spin, and therefore varies the degree of spin polarization 
of the electron liquid. At T = 0, for fields in excess of a critical field, H* ~ Ep, the electron gas is 
fully polarized, and consequently Ep is increased to twice its zero field value. The low temperature 
magneto-resistance in this case comes from a combination of the -fTy dependencies of the Fermi 
velocity, the density of states and the electron-impurity cross-section. In a FL, the doubling of Ep 
typically should lead to a decrease in the resistivity because scattering cross-sections are generally 
decreasing functions of energy. 

d. //|| dependence of the slope K = dp{H\\) / dT : At large magnetic field ( i/y > H*), the slope 
dp{T)/dT is significantly smaller than when ify = - by as much as a factor 100! This effect 
. las been observed in Si MOSFET's (See Fig. 3 in Ref. Jill) and in p-GaAs (See Fig. 2 in Ref. 

islp. In the FL framework there is no reasons for such a dramatic effect. Any scattering process 
responsible for the T dependence of the resistivity at H\\ = would typically be expected to give 
rise to T dependent scattering even when the electron gas is fully polarized. 

e. "Metallic" T dependence in samples with large ("insulating") resistivity, p > h/e^: The 
conventional theory of localization in the strong disorder limit p ^ h/e^ predicts that the electronic 
transport should be due to hopping conductivity. Therefore the resistivity of such conductors 
should increase with decreasing temperature and diverge as T — > 0. On the other hand, in some 
p-GaAs samples (See for example Fig.l in Ref. ^^) with large r^, the low temperature resistivity 
exhibits a "metallic" temperature dependence, i.e. p increases with increasing temperature, even 
when p > h/e^. 

f. Oscillations of p{H±) at T > Ep . Oscillations of the resistivity, p{H±) have been observed 
as a function of the strength of a perpendicular magnetic field, H± in both the Si MOSFET's and in 
p-GaAs hetero junctions. At relatively small r^, these oscillations were interpreted as Shubnikov-de 
Haas oscillations. However, for samples with large reported in Ref. jl^], these oscillations persist 
up to temperatures which are significantly larger than the bare Fermi energy. It is important to 
stress that FL theory not only predicts the existence of magnetic oscillations with period inversely 
proportional to the area enclosed by the Fermi surface, it also predicts that these oscillations should 
be entirely quenched when 27rT exceeds the cyclotron energy, ujc = eB/mc, which in turn must be 
less than or equal to Ep. Thus, the persistence of the oscillations to such high temperatures is in 
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disagreement with Fermi liquid theory ! 



B. Experiments on double-layers 

Additional information concerning correlation effects can be obtained from measurements of 
the "drag" resistance in a system of two 2DEG layers which are electrically unconnected. Current 
/ is passed through the lower (active) layer and the voltage Vd is measured on the upper (passive) 
layer. The drag resistance is defined to be the ratio pD = Vd/I- 

For relative small r^, experiments 0] on the drag resistance, pd{T), in double layer 2DEGs 
are in qualitative agreement with Fermi liquid theory [4^ . . Specifically, the drag resistance is 
small in proportion to (T/Ep)'^ and in proportion to {kpd)~°''^ where kp is the Fermi momentum, 
d is the spacing between the two layers, and typically = 2 or 4, depending on the ratio of 
d/L However, experiments [t?! on p-GaAs double layers with ~ 20 — 30, yield results which 
differ significantly from the predictions of Fermi liquid theory. These experiments are performed 
on samples with small resistances p ~ 0.05 — 0.1 h/e^^ in which quantum interference corrections 
to the Drude conductivity are insignificant: 

a. The magnitude of the drag resistance: The drag resistance in these samples is 1-2 order of 
magnitude larger than expected on the basis of Fermi liquid theory. 

b. T dependence of pd(T): Whereas in a Fermi liquid, Pd(T) ~ T^, in large devices 
Pd ~ (T)"-'^ where the temperature exponent exhibits non-Fermi liquid values 2 < ax < 3. (For 
example, ot = 2.7 in Ref. [131 •) 

c. dependence of pd{T,H^^): At low temperature, po{T,H\\) increases as a function of -ffy 
by a factor of 10-20 and saturates when H\\ > H* (See Fig. 1 in Ref. I3] )• In contrast, in the 
framework of the Fermi liquid theory, one would instead expect a significant decrease of pd{H\\) in 
the spin polarized state ( > H*) for the familiar reason that Ep is increased by a factor of 2 at 
H\\ > H* , and hence the electron-electron scattering rate decreases. 

A parallel magnetic field also appears to strongly suppress the temperature dependence of 
Pd{T) jl3| . In the presence of non-zero H^^, the value of aT{H\\) decreases with //y and saturates 
for H > H* at a value which is significantly smaller than the FL value ay = 2. (For example, in 
Ref. 17'] aT{H > H*) ~ 1.2.) 

d. The relation between p{T^H\\) and p£,[T,H\\): The T and especially the i/y dependencies 
of pd{T, H\\) and the resistivities of the individual layers p{HpT) look qualitatively similar to one 
another, which suggests that both have a common origin. (See Figs.l a,b in [l7| ). In a Fermi 
liquid, p{T, H^^) is primarily determined by the electron-impurity scattering, and Pd(T, //y) by the 
inter-layer electron-electron scattering, so there is no reason for their T and i^y dependencies to 
be similar. 
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C. Comparison with small Ts devices 

The anomalies discussed above have been observed in samples with large r^. In relatively smaller 
high mobility devices ( i.e. Vg ~ 1), behavior much more in line with the expectations of FL 
theory, modified by weak interference effects (See, for example, Ref. 37 1 for a review and Ref. [lo| 



for experiments in Si MOSFET's with high electron mobility and large distance to the gate. ) 

II. THERMODYNAMICS OF ELECTRONIC MICRO-EMULSIONS 

The two well understood phases of the 2DEG are a kinetic energy dominated FL phase (r^ ^ 1) 
and an insulating WC phase (r^ 1), which is a state of spontaneously broken translational 
symmetry, [o^l Were there a direct transition from a WC to an isotropic liquid, there are compelling 



reasons to think it would be first order. According to Landau, (See for example 84].) liquid 



crystal transitions are expected to be first order because of the presence of cubic invariants in the 



free energy functional. Brazovskii 33] showed that, even in cases (such as the liquid to smectic 
transition) where there are no cubic invariants, a freezing transition is expected to be fluctuation 
driven first order. However, we have proven (See below.) that, in general, no first order transition 
is possible when there are long, or moderately long-range interactions: V{r) ^ Q/r^ with D —\ < 
X < D + \ where D is the number of spatial dimensions. This result applies, therefore, to both the 
Coulomb (x = 1) and dipolar (x=3) cases in D=2. 

The apparent contradiction - that the transition must be first order and cannot be - is resolved 
by the existence of intermediate phases; there is no direct transition between a Wigner crystal 
and a uniform fluid phase. Rather, electronic microemulsions (ME), which can be visualized as 
consisting of intermingling regions of the two extremal phases, occur at intermediate densities. In 
essence, these phases are the vestige of the two-phase coexistence that would have accompanied 
the putative first order transition. The fraction of the system that is locally crystalline or fluid, 
as well as the size, shape and ordering of these regions are all determined by minimizing the free 
energy. It is important to stress that, since the energy associated with the patterns is extensive, 
these microemulsions constitute new phases in the thermodynamic sense. 3^ They provide an 
important paradigm of new non-Fermi liquid phases. On the mean field level, ME phases have 
many features in common with domain structures in 2D uniaxial ferromagnets [sil . Is^ and lipid 
films on the surface of water [5^. They are also similar to the various micro-phase separated states 
found in transition metal oxides 

0, H 0, 0, 0, 0, 0, 0, 0, S 

(including the stripe 



phases in the cuprate high temperature superconductors), the Bechaard salts|6ll. l62l. l63l|. and the 
ET based organic superconductors [g^. 



The mean field phase diagram for microemulsion phases has been discussed in Refs. |l9l.l2Cl.l21 1. 
This section therefore mostly describes the results, with the derivations (which generalize our earlier 
work) relegated to Appendix 1X1 The central result in the present section is Eq. |^ which gives the 
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FIG. 1: The T = phase diagram of the 2DEG in an "ideal" MOSFET, where d is the distance to the 
ground-plane. nc{d) is the putative critical density at which the free energies of the uniform Wigner crystal 
(WC) and Fermi liquid (FL) phases would be equal in the absence of microemusion (ME) phases. The solid 
lines mark the boundaries of the regime of the intermediate ME (stripe or bubble) phases. At mean field 
level, these lines are Lifshitz transitions. The hatched area represents the regime in which the regions of the 
two coexisting phases have sizes of order the electron spacing, so quantum fluctuations are order 1. ac is 
the critical value of d/as as discussed in Sec. IIIBI 

free energy density in the presence of either Coulomb or dipolar interactions of a microemulsion 
phase with domains of arbitrary shape and size. The mean-field phase diagram is derived by 
minimizing this equation. The resulting T = mean field phase diagram for the 2DEG in the 
dipolar case {e.g. a disorder-free MOSFET) is shown in Figs. ^andlS) The phase diagram for the 
2DEG without a ground-plane is obtained from this by taking the limit d — > oo, where d is the 
distance to the ground-plane. 

At finite T, the Pomeranchuk effect imphes that the fraction of WC increases with increasing 
temperature and //y , as discussed in Subsection 111 CI The resulting phase diagram is shown in 
Figs, m and 01 



A putative first order transition from the WC to the FL would occur when the free energy 
densities of the two uniform phases, FwciiT-) ^-nd FFLin-), cross at a critical density, Uc, as in Fig. 
m Recall that in the case of sufficiently short-range interactions, where such a transition can truly 
happen, it proceeds somewhat differently if considered as a function of chemical potential, fi, or of 
the mean density, n: As fi is varied through the critical value, the system jumps from a uniform 



A. Mean-Field theory of Electronic Microemulsions 
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FIG. 2: Mean field phase diagram of tlie 2DEG in an ideal MOSFET in the T - n plane. The dashed 
line indicates the putative phase transition between uniform FL and WC phases. The dotted line indicates 
the n-dependence of the bare Fermi energy. ME stands for microemulsion phases, rid is the lower critical 
density (proportional to d^^ where d is the spacing to the ground plane) at which the putative reentrant 
WC to FL transition occurs. 



WC with n = n_, to a uniform FL with n = n+, where n± are determined from Fwcin) and 
Fpi[n) by the Maxwell construction. Alternatively, as a function of n, the equilibrium state with 
n_ < n < n+ is macroscopically phase separated into two regions with density n_ and n+ and 
areal fractions {rij^ — n)/{n^ — n_) and (n — ?7-_)/(n+ — n„), respectively. 

If the interactions are of intermediate range, 2 < x < 3 the macroscopically phase separated 
state still has lower free energy than the uniform density states. (For x < 2, any macroscopically 
inhomogeneous state has infinite free energy density.) 

For present purposes, the reason interactions with a range x < 3 differ fundamentally from 
shorter-range interactions is that the surface tension is nonlocal, and is negative if the characteristic 
size of the domains, L, is big enough. In other words, one can still express the free energy of a 
general microemulsion in terms of a double integral over the interfaces between domains, but the 
integral is non-local. Specifically, any pattern of domains can be represented by a set of directed 
non-intersecting loops corresponding to the boundaries between the two phases, with the direction 
assigned so that when walking along the boundary, the WC phase is always to your left. In 
Appendix by a slight generalization of the calculations in Ref. 0, we show that the resulting 
free energy density can then be written as 

Fme = i2[fwc - fPL] + J dlao, -Jdl- dl'V{\l - I'D (1) 

where fwc = 1 — fpL is the areal fraction of WC, the integrals over dl are along the interfaces, 




FIG. 3: Mean field phase diagram of the 2DEG in an ideal MOSFET in the H\\ - n plane. The dashed 
line indicates the putative phase transition between uniform the FL and WC phases. Symbols are is in Fig. 

m 

and V{r) ~ Q^^r^'^^^'L (An explicit expression for V in terms of V is given in Appendix 1X1) 
Here do = <Jq{9) is the contribution to the surface tension from the short-range parts of the free 
energy, where 6 is the orientation of the interface relative to the a symmetry axes of the adjacent 
WC phase. Were V{r) shorter range than 1/r {e.g. for x > 3), there would be a well defined local 
surface tension, a'^-l'f = ctq — so whether or not there are microemulsion phases would 

depend whether or not cj'^-^-^ is negative. However, for the Coulomb and dipolar cases (and all cases 
in between), the surface tension is scale dependent and negative on large distances. Consequently, 
the free energy per unit length associated with long enough segments of interface is always negative, 
so the occurrence of microemulsion phases is universal. 

It is striking that V has the same power law behavior for x = 2ziz6x, and in particular for the two 
physically important cases of x = 1 and x = 3. The state which competes with the microemulsions 
is the uniform state for 1 < x < 2 and the macroscopically phase separated state for 2 < x, but 
the mean-field theory which determines the size and shape of the regions of the two phases that 
make up the ME is the same regardless of the sign of 6x. Specifically, in both the Coulomb and 
dipolar cases, Eq. ^can be expressed as 

Fme = - nc){fwc - Jfl) + | j dlao - j y|f=pp^^| 

where ai (for which an explicit expression is given in Appendix ^ is a microscopic coupling 
constant with units of a surface tension. We have included a short-distance cutoff, a, in Eq. |21 
which is roughly equal to the smaller of d and ^. The chemical-potential like term, /i, in Eqs. ^ 
and 121 is a function of the average density and the various interaction strengths. In the absence 
of the negative contribution to the surface tension from the long-range interactions, fwc would 
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FIG. 4: Mean-field n- dependences of tlie Free energy densities of the uniform FL and WC phases in the 
dipolar case. As is well known, in the Coulomb case the solid curves are convex at small n, but all that is 
important for present purposes is that the two curves cross with a finite relative slope dX n = tic- 

jump from its minimal value, fwc = 0, to its maximal value, fwc = 1 as /i passes from positive 
to negative. The negative effective surface tension is responsible for the existence of mixed phases. 
The mean-field phase diagram is found by minimizing Eq. ^with respect to the size and shape of 
the domains. 

Notice, however, that jl is not precisely a chemical potential, as it is conjugate to fwc rather 
than the mean density, n. In general, n as a function of /2 can be deduced by inverting the analysis 
that lead to Eq. [2| once the optimal solution of Eq. |2is found, the corresponding density profile 
can be then determined by means of Eq. IA14I in Appendix For the Coulomb case, the density 
profile in each region is non-uniform, and depends on the size and shape of the microemulsion. 
Specifically, the density deviates from its average value most strongly at the edges of the dipole- 
layer (of width ~ ^) associated with the WC - FL interface, and then decays toward the mean in 
proportion to the reciprocal distance from the interface. In the dipolar case, if L is sufficiently large, 
the mean-density within each region is essentially that of the corresponding uniform phase, i.e. 
n « f\Ycn+ + fpLn-. In this case, the values of fwc and fpL can be determined approximately 
by the Maxwell construction. 

One consequence of this is that, for the dipolar case is that /vyc ~ — — "--)) 

independent of the detailed nature of the ME, so the mean-field shape of the ME can be found 
by minimizin g iu st the surface free energy, second two terms in Eq. |21 at fixed fwc^ as was 
done in Refs. [20I . To obtain further insight into the meaning of the final non-local term, we 
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can compare it to the standard expression [80| for a finite size capacitor of area S and perimeter 
L:C = S/Aird + + 0{d). The same logarithmic term arises in integrating Eq. ^ 

over the boundary of an isolated domain, since indeed it has the same interpretation in terms of 
fringing fields. 



B. Electronic Microemulsions in an "Ideal" MOSFET 



Let us now consider explicitly the mean-field phase diagram of the 2DEG in the presence of a 
ground plane displaced by a distance, d, so that V{r) ~ e^/r for r <^ d and V{r) ~ e^d^/r^ for 
r ^ d. (The Coulomb case corresponds to the limit d ^ oo.) The mean- field phase diagram for 
this problem is shown schematically in Fig. ^ 

If there were only the uniform FL and WC phases, there would be a single line of putative first 
order phase transitions, indicated by the dashed line in Fig. ^ A peculiarity of this phase diagram 
is the re-entrant character of transition as a function of n. For Coulomb interactions, the ratio of 
the potential to kinetic energy, r^, decreases with increasing density, so where nc?^ ^ 1, increasing 
density always favors the FL. Conversely, for short-range interactions, as in ^He, is an increasing 
function of density, and hence increasing density favors the WC. From this point of view, dipolar 
interactions are much the same as any other short-range repulsions, so where nd'^ < 1, the FL is 
the low density phase. One implication of this analysis is that there exists a critical value of the 
ratio, d/as = etc, such that for d < OcOs, there is no WC phase at any density. In the same way 
that the value of at which the energy of the uniform WC and FL cross estimated in [s^ turned 
out to be large, we think that is likely that that Oc ^ 1, but as far as we know, this issue has not 
been addressed by any of the numerical methods that have been applied to the problem. 

The microemulsion phases exist between the two solid lines in Fig^ We will characterize them 
first at mean field level and later with fluctuation corrections included. At mean field level, we 
must minimize the free energy with respect to the size, shape, and density profile of the domains of 
the two phases. A complete analysis of the patterns that minimize Eg. El does not exist. However, 
the general structure of the solutions has been discussed in Refs. |l9|, |20|, |21 1 , as summarized in 
Fig. 

In the middle of the phase separated region, where fwc = I'fl, the free energy is always 
minimized by a stripe array. The integrals in Eq. |21 are easily done for stripes of width L 3> a, 
with the result 

Fstripe = L'^Wo - (Ji log(L/a)]. (3) 

This expression can be minimized with respect to L to obtain the equilibrium stripe width, 

L = Lq = aexp[ao/cri + 1]. (4) 

So long as Lq 3> a is large {i.e. if o"o > cri), the considerations that lead to this result are self- 
consistent. If Lq is comparable to the spacing between electrons, fluctuations and other effects 
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FIG. 5: A schematic picture of the mean field micro-emulsion phase diagram of an "ideal" MOSFET. The 
solid bars indicate regions of more complicated (super-micro-emulsion) phases, consisting of a mixture of 
stripe and bubble crystal regions. 

could qualitatively affect the nature of the equilibrium state. However, even in this case it still 
follows from the large L behavior of Eq. |21 that there is a ME state that has lower energy than 
the uniform or macroscopically phase separated states, even if the present mean- field theory is 
insufficiently microscopic to describe it. 

In a simple (but not generic) situation, where the surface tension is strongly anisotropic across 
entire phase diagram, there are only uniform and stripe phases |2o| . (This can readily happen, for 
instance, in semiconductor heterostructures grown along a cyrstalographic direction with a highly 
anisotropoic effective mass tensor.) In this case, as n tends towards n_, the stripes of FL continue 
to have a characteristic width of order Lq, while the intervening regions of WC have a diverging 
width as fwc ~^ 1- (The same story pertains to the limit n — > n^, with the role of the WC and 
FL regions interchanged.) We will refer to this type of transition (which recurs several times in 
our analysis) as a "Lifshitz transition," despite some differences with the usual Lifshitz case, jo^ 

When cr is less anisotropic, including the usual WC case when a has the same 6-fold rotational 
symmetry as the triangular lattice, the number of distinct phases in the mean-field phase diagram 
is much larger: In the limit when the fraction of minority phases is small, Fme is minimized by 
a bubble crystal phase. For n near n_|_, this phase consists of a triangular lattice of crystallites 
embedded in the liquid. Conversely, for n near n_, the bubble phase consists of puddles of FL 
embedded in WC. The distance between the bubbles diverges, but their radius approaches a con- 
stant, , Lb ^ Lq, as the fraction of the minority phase tends to zero, so again there is a "Lifshitz 
transition" to the uniform phase. 

It follows from Eq. |2that the interaction between bubbles at large separation R ^ Lq always 
is dipolar, as 



A direct transition from a bubble to a stripe phase at a critical value of fwc would, at least at 
mean-field level, be first order, and hence is forbidden. Iterating the above analysis, this putative 
transition must give way to a new set of super-microemulsitons, consisting of regions of bubble 
phase and regions of stripe phase. These phases are indicated in Fig. El by black regions. Since the 




(5) 
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surface tension of a stripe phase is highly anisotropic, it is probable that one super-stripe phase 
(as opposed to an infinite hierarchy of them) result from this analysis, as has been discussed in 
Ref. y 

The various possible ME phases can be grossly classified by their broken symmetries. They can 
also be classified by whether they are insulating or fluid. Man y o f the phases can productively 
be thought of as quantum analogues of classical liquid crystals [70|: For instance, a stripe phase 
which breaks translational symmetry in one direction, only, but conducts readily in the transverse 
direction is analogous to a smectic liquid crystal, and therefore can be referred to as an "electron 
smectic." This analogy also permits us to conceive of phases that could be thought of as partially 
melted versions of the various mean-field phases: A partially melted stripe phase, for instance, can 
naturally give rise to an "electron nematic" phase [toI . Iq^ . 93 1 in which translation symmetry is 



restored, but the preferred stripe direction remains distinct, leading to a uniform fluid state with a 
spontaneously chosen nematic axis. For the most part, we will leave the exploration of these many 
rich possibilities to a future study. 

C. Thermal physics: the Pomeranchuk effect 

The thermal evolution of the microemulsion phases at low temperatures is dominated by the 
spin degrees of freedom [19I . I20I 1 . In direct analogy with the Pomeranchuk effect familiar from the 
theory of ^He, the spin entropy of the WC is large compared to that of the FL, so the system tends 
to freeze upon heating! 

In the FL state, the electron spin degrees of freedom are quenched due to the Fermi statistics. 
Consequently, the Fermi liquid entropy density is small, Spi ~ ^{T/Ef), at T <C Ep. Conversely, 
at T > J the spin entropy of the WC is large. The exchange interaction between localized electrons 
is exponentially small, J ~ exp(—Qy^). (See, for a review, 0]). For example, a WKB calculation 
for Tg = 30 yields J ~ 5mK. Therefore, at essentially all experimentally accessible temperatures 
(as of this date), T > J. Consequently Swc ~ nln2 ^ Sfl, and 

FwciT) ~ FwciO) - Tn log [2]. (6) 

This results in a T dependent shift in the "center" of the microemulsion phase, nc{T): 

ne(T) =n,(0)[l + Cr'riog(2)], (7) 



where Ci is defined in Eq. IA21 As a result, anywhere in the microemulsion phase, fwc increases lin- 
early with T. This behavior will dominate the T dependence of many important physical properties 
of the system at low T. 

At T > Ep the liquid is not degenerate and therefore there is none of the entropy suppres- 
sion characteristic of the FL. Only at these high temperatures does the liquid entropy exceed the 
WC entropy. One consequence of this is that for densities n ~ ric, the WC melting temperature 
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Tm ~ Ep. Thus there is a single temperature scale, T ~ Ep, that governs the crossover to a low 
temperature electronic microemulsion regime. The effect of temperature on the phase diagram of 
electronic microemulsions is shown qualitatively in Fig. [21 Of course, within the microemulsion 
regime there may be a number of thermal phase transitions between states with different pat- 
terns of symmetry breaking, such as stripes, bubbles, nematics, etc. These more detailed thermal 
transitions will be addressed in the next susbsection. 

Since the Pomeranchuk effect is the physics of spin entropy, similar considerations govern the 
-ff|| dependence of the phase diagram (See Fig. A parallel magnetic field couples only to 

the electron spins. Since the spin susceptibility xwc ^ Xfl, the corresponding magnetization, 
^wc ^ at small H^^. Since the free energy of the system contains the term —MH^^ there is 
a -fTy induced increase of the WC fraction over a wide range of circumstances. 

For Tm ^ T > J, the combined //y and T dependence of Fwc is that of free spins: 

Fwc{T,H\\) = Fh/c(0,0) - TnS (8) 
5 = log[2cosh(/7||/r)]. 

In this same range of temperatures, the T dependence of Fpi is negligible, and so long as H\\ <^ Ep, 
its dependence can be neglected as well. In this regime, therefore, the fraction of WC can be 
expressed as a scaling function 

fwc{n, H\\ , T) = fwc{n, 0, 0) + gin, TS) (9) 
for T <^Ep k H\\4;:H* = Ep, 

where g is a, non- negative monotonically increasing (initially linear) function of TS, and g{n,0) = 
0. At high enough magnetic field, when all the spins are polarized, the Pomeranchuk effect is 
completely supressed, so that 

fwc{n,H\\,T) = fwc{n, 00,0) (10) 
for T <^Ep h H\\> H\ 

1. Experimental consequences of the Pomeranchuk effect 

Since the resistance of the WC and FL are very different, the Pomeranchuk effect is one of the 
most directly testable features of the microemulsion phases. The expected T and H\\ dependences 
of the resistance of a micrio-emulsion phase p are shown qualitatively in Figs. 7 and 8 in Sec. 
IVIIII The fact that p{T,H\\) increases with T and H\\ at low T reflects the fact that the fraction 
of WC increases as T and H\\ increase. The fact that at H\\ > H* , the T and H\\ dependences of 
p{H\\) are quenched reflects the fact that there is no spin-entropy in the fully polarized system, so 
fwc{T,H\\) no longer depends strongly on these quantities. At higher T, the WC melts and the 
decrease of p{T) reflects the corresponding T-dependence of the viscosity of the electron liquid. A 
more detailed analysis will be presented in Section. IIII CI 
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D. Fluctuation effects 



The mean field treatment gives a reasonable order description of the micro-emulsion phases 
when the characteristic spacial scale is large, toLq ^ 1, so that the elementary structures of the 
micro-emulsion contain many electrons. In this case, where the mean-field theory fails, it can only 
be due to the effects of collective long-wave-length fluctuations on length scales large compared to 
Lq, rather than due to large- amplitude fluctuations on length scales of order Lq. 

Let us begin by considering the case of the bubble crystal. Near the Lifshitz point, where the 
spacing between bubbles R ^ oo, the thermal energy of bubbles is larger than the interaction 
energy Eq. |S1 so the melting temperature of the bubble superlattice must vanish as i? — > oo. On 
the other hand, for fixed R, the strength of the interaction between bubbles is finite, so the the 
bubble crystal phase survives thermal fluctuations in the usual sense that the correlations of bubble 
positions exhibit power-law decay. The nature of the thermal transition between the bubble crystal 
phase and the uniform phase is not, presently, settled. Of course a direct first order phase transition 
is forbidden. One possibility is that there is a sequence of two transitions with an intermediate 
hexatic phase. Alternatively, there may be a further set of hierarchical microemulsion phases. 

Let us now discuss the role of the quantum fluctuations of the positions of the bubbles. At low 
bubble concentration, oc R'^, they can be characterized by an effective mass Mb- (See Sec. 
IV Al ) Thus we can introduce a ratio between the potential energy of the interbubble interaction 
and the zero-point kinetic energy r^-l^f ~ R^^Mb- Since r^^^ becomes small at small values of hb, 
at T = 0, the bubble crystal phase is unstable below a critical value of n^. Consequently, the 
Lifshitz point of the FL-WC bubble phase transition is preempted by quantum fluctuations. A 
quantum melted WC bubble phase has the same symmetry as a liquid. The nature of the quantum 
phase transition between the bubble liquid and bubble crystalline phases is currently unknown. 

Let us now consider the role of thermal fluctuations on the stripe phases. If the Hamiltonian 
is rotationally invariant then the smectic phase is unstable at any non-zero temperature to the 
proliferation of dislocations. Thus the mean field smectic phase is replaced by a nematic phasejoJl) 
which in keeping with the Mermin- Wagner theorem, does not actually break rotational symmetry. 



but rather has power-law orientational order. We refer the reader to Ref. 20| for details. 



E. Crossovers at higher temperatures 

Although it is beyond the primary focus of the present paper, it is worth commenting on the 
nature of some of the higher temperature behavior of the electron fluid at large r^. 

At very low electron densities and T = 0, the WC is classical, and there are two characteristic 
energies in the system: the electron interaction energy V = e^-y/vm and the plasma frequency 

Saluated at wave- vectors of order n^/^), Op ~ \J e^np/'^ /m <^ V. The classical melting occurs 
, 3| at a critical temperature T^i ~ V/125 where V ^ T^i ^ Qp. For T > T^u there are still 
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FIG. 6: The phase diagram of the 2D electron gas. FL, WC and ME stand for the Fermi hquid, Wigner 
Crystal and electronic micro-emulsion phases respectively. The dashed lines indicates the n-dependencies of 
the Fermi energy, the plasma frequency, and the inter-electron interaction energy. 



two distinct temperature intervals: 

1. V > T > Tci: Here the electron fluid is a highly correlated classical fluid. There are many 
examples of such fluids - indeed, most fluids fall in this regime js^l- In our opinion, the present 
theoretical understanding of such strongly correlated fluids is rudimentary, but there is a wealth of 
empirical information that we can draw on to make fairly confident statements about the expected 
behaviors. In particular, the viscosity is widely observed 50, 6^ to be an exponentially increasing 



function of decreasing temperature, r] ~ exp[ro/T] where Tq ~ V . 

2. T > V: Here the electron system forms a classical weakly interacting gas, where the viscosity 
is an increasing function of T because the electron mean free path is an increasing function of energy 

Q 

Let us now consider the case when the electron density is near the critical density of the putative 
quantum FL-WC transition. The critical value of = r* is "of order one". However, it turns 



out 



53] to be numerically quite large: r* ~ 40, so we hope to be forgiven for considering it as 



a large parameter. In this case there are three characteristic energies: 1) Ep, or more properly 
25l . |26| . the renormalized Fermi energy, Ep < Ep, which contains a renormalized mass, 2) the 
interaction energy V = rgEp, and 3) the plasma frequency ilp ~ \fEpV = y/rlEp- For large r^, 
these energies are quite distinct, with V > VLp > Ep, as shown schematically in Fig. El Thus, 
we can identify three distinct regimes of temperature, above the transition temperature to the 
microemulsion phases, where the electron fluid can be expected to exhibit different qualitative 
behaviors: 

1) T < Ep-. Here, of course, as long as the system is fluid, quantum coherence is important as 
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are the Fermi statistics. 

2) Qp > T > Ep-. Here the electron-fluid is in a semiquantum regime |25l. . Quantum coher- 
ence, and specifically the Fermi statistics of the electrons, are presumably not terribly important 
[ssl ]. The motion of electrons consists mostly of quantum oscillations around certain positions, 
with a frequency fip ^ T, which is also large compared to the rate, l/rcage, at which electrons 
escape from the local "cage" formed by the neighboring electrons, fip S> 1/Tcage- As far as we 
know, there is no theory of the behavior of electronic fluids in this reg ime. However, a similar 



regime has been discussed in the framework of liquid "^He and '^He 

that l/Tcage ~ Ep. 

3) V > T > Qp: Here the electron fluid is a highly correlated classical fluid. 
4:) T > V: Here the electron fluid forms a classical gas. 



26l |. where it was argued 



III. BEYOND BOLTZMAN TRANSPORT IN CLEAN, CORRELATED ELECTRON 

FLUIDS 

This section is something of a digression. The issue addressed here is how the resistivity of a 
highly correlated fluid in the presence of weak disorder is related to the hydrodynamic properties 
of the ideal fluid in the absence of disorder. The considerations here are general, and apply to 
a FL as well as to an electronic microemulsion. Despite this, much of the analysis has not been 
carried out previously, to the best of our knowledge. So we will discuss the general problem in 
the present section, and then apply the results to compute the transport properties of electronic 
microemulsions in the following sections. 

The resistivity, p, of the electron gas is determined by the electron-phonon, electron-electron 
(ee) and electron-impurity (e-i) interactions. At the low temperatures of interest, electron-phonon 
scattering is negligible. Since kp is much smaller than the reciprocal lattice vector of the atomic 
lattice, e-e Umklapp scattering is also negligible which means that e-e scattering processes con- 
serve quasi-momentum and so do not directly contribute to the resistance. Thus, ultimately, the 
resistivity is determined by the electron- impurity interactions. However, the actual mechanism of 
momentum transfer from the electron system to the lattice can vary depending on circumstances. 
In the majority of cases considered in the literature (see for example [29|) the current distribution 
is treated as statistically uniform and the e-e and e-i scattering are treated as uncorrelated, and 
are described by different parts of the scattering integral in the Boltzmann kinetic equation. (See 
Subsection IIII Al ) In this case, e-e scattering can contribute to the resistance only to the extent 
that it affects the form of the electron distribution function. In particular e-e scattering does not 
contribute to the resistance at all if the electron-impurity scattering rate is energy independent. 

In Subsection IIIIBl we show that the assumption concerning the uncorrelated nature of e-e 
and e-i scattering is not valid when the electron electron interaction is strong and the disorder 
is sufficiently weak. The main result of this subsection is that in this limit, the resistance of the 
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sample is proportional to the electron viscosity, 

poc-q. (11) 

To appreciate how significantly this result differs from the usual Drude result, note that rj in a FL 
is proportional to the e-e mean-free path, i^e- (See Subsection IIII CI ) Since iee is a decreasing 
function of T, this implies the opposite sign of the temperature dependence of typical metals! More 
fundamentally, the strong interference between correlation effects and impurity scattering implied 
by this relation is impossible in the context of the Boltzman equation. 

We derive Eq. ^Jin several limiting cases in the following. The simplest case (See Subsection 
HUB 1\ ) is where there is a dilute set of large radius obstacles placed in an otherwise ideal FL. Then 
the force on the obstacles from the moving liquid is given by the Stokes formula and is proportional 
to r]. The important point is that, under any circumstance in which the disorder leads to a non- 
uniform distribution of the current on scales long compared to the appropriate microscopic healing 
length (iee in the case of the FL), the resistance is more or less proportional to rj, as in Eq. 1111 and 
so will exhibit T and H^^ dependencies that reflect the corresponding dependencies of the viscosity 
r/(T, ) in the ideal system. 



A. Boltzman-Drude Resistivity 



Usually electron-electron and electron-inpurity scattering in semiconductors and metals are 
described by a Boltzmann kinetic equation 

E • Vrnk(r) = -Iee(k, [nk]) - hiiK M (12) 

where /ge and lei are electron-electron and electron-impurity scattering integrals which are func- 
tional of the electron distribution, nk(r), but are independent of r; and E is the electric field. 
One of the assumptions which leads to Eq. is the absence of correlations between the e-e and 
e-i scattering, which in turn is related to an assumption that in a statistically uniform system the 
distribution function of the electrons nk is spatially uniform, i.e. it is the same near impurities 
and between them. 

Since the e-e scattering conserves the total quasi-momentum of the electrons, only electron- 
impurity scattering contributes to the resistivity leading to the Drude formula (See for example 

m 

PDR = — (13) 

where v is the density of states at the Fermi level, vp \& the Fermi velocity, I/tj is the electron- 
impurity transport scattering rate computed in the usual way. 

If the impurity potential is treated in the white noise approximation, Tj is independent of the 
electron Fermi wave length. Thus p{T,H\\) is independent of T and H\\. Typically, though, the 
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electron-impurity cross-section is a decreasing function of the energy so p{T,H\\) decreases both 
with T and with H\\. 

As is weh known, the Drude-Boltzman description breaks down at low enough T, due to quantum 
interference effects ; "weak localization" phenomena is an example of this sort of corrections. 
However, in the weak disorder limit, these corrections are only of large magnitude at exponentially 
low temperatures. 

What we discuss below is a different sort of breakdown of the Boltzman description, which is 
important in the limit of small disorder, especially in the case of low-dimensional highly correlated 
systems. Specifically, the assumption of Boltzman theory that the current distribution is uniform 
is valid only so long as the impurities responsible for the resistivity are sufficiently dense that they 
can be treated as self-averaging. Roughly, this requires that electron-impurity mean free path is 
smaller than electron-electron one iei ^ee- 



B. Relation between resistance and viscosity 



For sufficiently weak disorder, the motion of a FL at length scales larger than 
in the context of the hydrodynamic Navier-Stokes equations [35 1 



n 
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dt 



+ Ui 



dUi dF 



drj ~^ dvi 
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dri 



can be treated 



(14) 



where Ui{r) are the components of the electron hydrodynamic velocity, F is the free energy density, 
Iljj is the momentum current tensor, and ^*(r) is an external potential (i.e. the impurity potential). 
Hij, in turn, can be expressed as 
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6 dri J dri 



(15) 



Here are the first and the second viscosities respectively. The first viscosity r]i determines the 
viscous drag on an object moving through the fluid, and so is the more important one for present 
purposes. The second viscosity 772 is associated with the rate of equilibration of density variations. 



1. Dilute large obstacles 



The first case we will consider is one in which the "disorder" consists of a dilute concentration, 
Ni of large obstacles (radius A ^ l^e but still NiX^ <^ 1) in the electron liquid. (For example, 
these obstacles could be pinned bubbles of WC, as we will discuss in Subsection IIVB . below.) The 
resistance can be computed straightforwardly by equating the electric field to the force exerted on 
a flowing liquid by a stationary object, which is given by the Stokes formula. As a result we have 
an expression for the resistivity 

P = , ^ 1/2., ' (16) 

|ln(AiV//2)| 
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The fact that p depends only logarithmically on A and is a slightly superlinear function of Ni 
reflects the "Stokes paradox" in 2D hydrodynamics 

2. Dilute microscopic impurities 

The case of a 2D FL with i^e larger than the electron- impurity cross-section, but still small 
compared to the spacing between impurities was analyzed in Ref. 3J]. The result is: 

PDR 



[l + ^ln(/eeiV//')] 



(17) 



where a is a number of order 1. The logarithm in the denominator of this formula is also a 
consequence of the Stokes paradox in 2D. Note that when the logarithmic term is large, this again 
gives a resistivity that is proportional to the viscosity. 

3. Weak, smooth disorder potential 

Another case which can be analyzed readily is that in which the scattering potential is smooth 
and is characterized a length scale of A ^ Zee and of small amplitude Vq . 

(18) 

C. Hydrodynamics of an "ideal" electron liquid 

Before embarking (in the following sections) on the analysis for the micro-emulsion phases, in 
this section we review briefly what is known about the viscosity of FL's. We also discuss the vis- 
cosity of strongly correlated electron liquids above the degeneracy temperature. Surprisingly, there 
are still many regimes in which the answers are not known, either theoretically or experimentally. 
We therefore also invoke experimental results on model liquids, especially ^He, and in one case 
even ^He, to support a phenomenological conjecture concerning the T dependence of the viscosity. 

1. Viscosity of a pure Fermi liquid 

In the Fermi liquid regime (T ^ Ep) the two viscosities are of the same order (?/fl ~ ~ ^2 
and 

riFL ~ npFVpTee (19) 



is inversely proportional to the square of the temperature and becomes very big at small T 29|, 
because the electron electron scattering rate in the Fermi liquid 

T-}^C— for T-€.Ef (20) 
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decreases with T. Here C = C^^ + C|j is a coefficient of order one which is proportional to the 
sum of scattering probabihties for parallel and anti-parallel spins. 

To characterize the -ffy dependence of the viscosity we discuss the ratio 

r?FL(0) ^ (ree/m*)j^,|=o ^ (C||m*)j^|| 

r]FL{H\\) (2Tee/m*)H|| ^ 4(Cm*)H||=o ^ ^ 

evaluated for //y > At high electron concentration when the interactions can be treated 

peturbatively (r^ <C 1) the value (r/(0)/r/(f/'|| > H*) < 1: both the scattering amplitude and the 
electron effective mass decrease with increasing spin polarization of the electron system, and the 
Fermi energy increases. In the case > 1 the situation is more complicated and currently it is 
unknown whether m*{H^^) increases or decreases with H^^; the Fermi energy, however, clearly still 
increases with H^^. A similar situation (large Fermi liquid) pertains to ^He, but the dependence 
of rjFL on spin polarization is unknown there, as well. 



2. Viscosity of a nondegenerate strongly correlated electron liquid 

Eqs. El and 1221 are valid in a FL at any as long as the temperature is low enough and the 
system is in a liquid state. However, if > 1, then (as mentioned above) there is a temperature 



interval Tp <^ T <^ Qp = yJEpV where the fluid is still quantum but not degenerate, and is 
still strongly correlated. The temperature dependence of the viscosity of the electronic liquid in 
this interval is unknown. However, a similar regime exists in another strongly correlated quantum 
liquid: ^He. We know of only one theoretical paper 
presented that 

1 



26l | on the subject where the arguments were 



r/ ~ - for Ep <T < ^EpV . (22) 



Since the theoretical arguments made in Ref. [2a| are quite general it is plausible that the T- 
dependence of the electronic viscosity in the semiquantum regime can be described by the same Eq. 
1221 . Surprisingly, even experimental data to test this conjecture in ^He is unavailable. In support 
of it, the author of Ref. ^23] observed that the viscosity of '^He in the same temperature range 
is roughly consistent with a 1/T dependence; since in this range of temperatures the statistics 
of the particles is likely unimportant, this evidence should be taken as moderately compelling. 
This means that, as in classical liquids, the viscosity of the semiquantum liquid decreases with 
increasing T. The difference, however, is that it decreases only inversely in proportion to the 
temperature, whereas the viscosity of classical correlated liquids decreases exponentially with T 



5C 



60|. Presumably, the same is true for the electron- fluid 



7/~ exp[-Ay/T] for ^/EpV <T <V , (23) 

where ^ is a number of order 1. 

Lastly, for T > y we have is a classical electron gas where the viscosity increases with T. 
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IV. THE NON-DEGENERATE BUBBLE LIQUID 



The liquid phase consisting of a fluid of WC bubbles in a sea of FL is the simplest of the 
microemulsion fluid phases. We will start with the simplest sub-case in which the system is near 
the mean field Lifshitz point, where the bubbles are weakly interacting, and where (as we showed 
in Section III D|) even small thermal and/or quantum fluctuations destroy the bubble-crystalline 
order. We moreover confine ourselves to the limit in which the bubble-size is large compared to the 
spacing between electrons, nL'^ ^> 1, and consequently the bubbles are well defined as droplets of 
WC, and are relatively heavy so that quantum coherent motion of the bubbles occurs only at very 
low T. Slow and long wave-length motion of the electron liquid can be described with the help of 
the Navier-Stokes Eqs. 1141 and IT^ However, the viscosity of the system is very different from the 
viscosity of the FL. In particular, the viscosity depends on the bubble concentration, ub- Thus 
(following Eq. the resistivity in the presence of small disorder directly reflects the T and i^y 
dependencies of fwc discussed in Sec. ^ above. 

Because rj depends on T directly, and implicitly through the T dependence of ub, the final 
result is complex and there are many special cases, as we shall see. One important result is that, 
in contrast to the FL case, the viscosity extrapolates to a finite value at T ^ 0. Often, is a 
non-monotonic function of T. Near the boundary of the bubble liquid phase, where is large, rj 
typically depends exponentially on T. 



A number of scattering processes contribute to the hydrodynamics. Within the majority FL 
portion, the same electron-electron scattering process, r^e, that determines the viscosity of the 
uniform FL plays a major role in the bubble liquid, as well. However, there is also an important 
effect of the scattering of electrons (quasi-particles) of the FL off the WC bubbles, which occurs 
with rate 



where Lb Lq (defined in Eq. |^ is the bubble diameter. Although the bubbles are dilute, they 
are strong scatterers. 



The effective scattering rate which determines the viscosity is the sum of contributions of dif- 
ferent processes, 



A. The hydrodynamics of the ideal bubble liquid phase. 



1/TeB ~ VF TIbLb 



(24) 



1. 



First viscosity of the dilute bubble liquid: 



1 




(25) 
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At low temperatures r^e ^ Te_B, the electron-bubble scattering is the dominant scattering process. 
Consequently, 

rii ~ nkp — (26) 

i.e. the first viscosity approaches a T independent constant. This is very different from Fermi liquid 
temperature dependence Eq. according to which rj diverges at T ^ 0. 

At higher temperatures, where fee is less than the distance between bubbles, the viscosity is more 
or less dominated by the electron-electron scattering, as implied by Eq. 1251 although with some 
quantitative corrections due to the presence of bubbles. For instance, at high enough temperatures 
that /ee *S Lb, the hydrodynamic description works even on distances smaller than the inter-bubble 
distance. In this case the viscosity follows the usual formula for viscosity of classical suspensions 

r]i = r,FL{l + nBLl) (27) 

(still assuming that ubL^ ^1). 

The above analysis ceases to be valid at very low temperatures where the motion of the gas 
of bubbles itself becomes quantum coherent, as discussed in Section below. Nevertheless, for 
sufficiently dilute bubbles, the crossover to a quantum coherent bubble fluid occurs at such low 
temperatures that it is not worth belaboring. 



2. First viscosity of a dense bubble liquid 



Let us consider now the interval of temperatures and densities where the system of bubbles is 
reasonably dense and close to forming a bubble crystal. The situation here is similar to that in 
the theory of classical liquids: the complete theory is not available, but it is plausable that the 
flow is associated with local rearrangements of the bubble configurations. Thus, in this regime, as 
discussed in Section HTll 

m ~ exp[VB/T] (28) 
where Vb is a characteristic bubble-bubble interaction energy. 



3. Magnetic field dependence of the first viscosity 

At low temperatures, when Eqs. UHl and |27| are valid, the i?||-dependence of the viscosity is 
determined primarily by the corresponding H^^ dependence of ub- Since fwc{H\\) = nB{H\\)L\, 
and Lb is only weakly iify dependent, fwc increases linearly with H^^ at small //y, and saturates 
for Hii > H\ 
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In the case of large bubble concentrations, the bubbles form a strongly correlated suspension 
and Eq. |2H1 holds. Since ns increases with increasing presumably Eq does as well, so the 
viscosity is an exponentially increasing function of i^y . 

4- Second viscosity of the bubble liquid 

It is well known that in the presence of a loiig internal relaxation time in the liquid, the second 
viscosity can be much larger than the first one [35|. In the case of the WC bubble phase the origin 
of a long relaxation time is the following: As the total electron density changes, the equilibrium 
bubble concentration, ub, changes as well. On the other hand, the nucleation of a new bubble 
occurs at an extremely small rate, with an activation energy (for thermal nucleation) or a tunnel 
barrier (for quantum nucleation) which is grows with increasing Lb- For example, the height of 
the thermal nucleation barrier can be estimated from Eq. ^ to be 

^0 ~ (^Lb, (29) 

leading to an exponentially long relaxation time, 

^2 ~ exp(^) (30) 

B. Resistivity of the WC bubble liquid with small disorder 

In the bubble phase there are two types of the current carriers: electrons and bubbles. The 
relaxational dynamics depends on the interactions of the electrons (quasiparticles) and bubbles 
with each other and with the impurities. Depending on the relative strength of these different 
interactions, there can be many distinct regimes of behavior. Here we consider only two. 

Case 1: Electron mobility 3> bubble mobility 

If the bubbles are big and massive, then under many circumstances their mobility will be much 
smaller than the mobility of the quasiparticles in the FL. In this far as the calculation 

of the resistivity is concerned, the bubbles can be treated as stationary objects which scatter the 
electrons. They differ from ordinary impurities in that their concentration varies as a function of 
T and H\\, and in that they can move. 

At low temperature, where £ee is long compared to the elastic mean free path, the resistance is 
given by the usual Drude formula /9 ~ 1/ Te/ / where 

^ = i + i. (31) 

Teff n TB 

and 1/tb oc ub- Thus the T and //y dependencies of the resistance are determined by the Pomer- 
anchuk effect, as in Eqs. |H1 We think that this is the case which is most realistically realized 
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in currently available experiments. The characteristic T and ify dependencies of the resistance 
are shown in Figs. 7, and 8 The saturation of the resistance at large -fTy is due to the complete 
polarization of the electron liquid. The downturn of the T dependence of the resistance at T > Ep 
takes place in the region where the Wigner crystal melts and the temperature dependence of the 
resistivity R(T) is determined by the temperature dependence of the viscosity, Eq. |2S1 

It is worth noting that the enhancement of the resistance by the scattering from bubbles is 
particularly dramatic at low temperatures in the case where the impurity potential is smooth. As 
is well known, a smooth potential causes little back-scattering. However, when the bubbles are 
large and have lower than average electron density, they will tend to get stuck on the mountaintops 
of the smooth potential. The bubbles are effective large angle scatterers. 

Case 2: Electron mobility comparable to bubble mobility 

In the opposite limit of clean samples and at finite temperature, the equilibration between the 
motion of the bubbles and the electrons takes place more quickly than the momentum transfer 
from the electron system to the impurities. In this case the resistance is determined by the Stokes 
hydrodynamic formula and the T and H^^ dependences are expressed in terms of the viscosity of 
the suspension, as in Section Hi CI 

V. QUANTUM BUBBLES. 

In this section we will discuss some of the issues concerning the role of the zero temperature 
quantum fluctuations of bubbles of a minority phase which are beyond the treatment of the previous 
section. Quantum motion of the bubbles is most important when the bubbles are small, nL^ ~ 1 
— the situation which applies, for example, in the shaded portion of the phase diagram in Fig. ^ 
However, we do not have the analytic tools available for treating this most interesting regime of 
parameters. 

When the bubbles are large, nL^ ^ 1, quantum fluctuations are small under most circum- 
stances, and consequently the ground-state is an ordered, bubble crystal. However, near the Lif- 
shitz point, where the density of bubbles is small, ubL^ <C 1, quantum fluctuations are sufficient 
to quantum melt the bubble crystal, even when the bubbles themselves are large. We therefore 
focus on this rather extreme regime as being a more tractable limit of the more general problem. 

In the Subsection IV Al we discuss estimates for the effective mass Mb of WC bubbles. In 
Subsection IV Bl we show that as long as the resistance of the sample is smaller than h/e^, the WC 
bubbles are not localized by frozen disorder inspite the fact that they can be heavy objects. 
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A. Bubbles as quasiparticles 

As discussed in Sec. IIID| at low density the WC bubbles form a quantum fluid. A WC bubble 
has a finite lifetime with respect to quantum melting, which, for a large bubble, is exponentially 
long. Thus, for large WC bubbles, there is a broad intermediate time interval in which they behave 
as propagating quasiparticles with a well defined effective mass Mb and a charge e* = eMB/m. 
(The latter estimate follows from the fact that in a liquid, the momentum and the flux of mass 
have the same form. ) 

There are two different estimates of the bubble effective mass one can make, depending on the 
quantum character of the interface between the WC and the FL: 

L%n quantum smooth surface , , 

f : (32) 

L^An quantum rough surface 

If the surface is quantum smooth, ( which means that in the ground state there is zero density 
of steps on the surface) then the bubble moves as a solid object and Mb is determined by a 
redistribution of the mass of the liquid over distances of order the radius Lb (upper formula in 
Eq. I32J. On the other hand, if in the ground state the surface has steps (it is quantum roug h Q) 
then the bulk of the crystal does not, actually, move. Rather, a strip near the rear quantum melts, 
and a stripe near the front coherently crystalizes. Only an amount of mass proportional to the 
difference in the densities of the WC and the FL is transported as the crystal translates by this 
process. Under this assumption, we arrive at the estimate given by the lower expression in Eq. 1321 

Now we would like to make few comments about the bubble liquid phase that occurs near the 
microemulsion to WC phase boundary. Here, the WC is the majority phase, and well ordered, but 
with dilute bubbles of FL inside it. FL bubbles embedded in the WC are topological objects, so 
the number of additional electrons associated with each is an integer. Thus the statistics can be 
either fermionic or bosonic depending on the parameters of the bubbles. 

We would also like to mention the role of electron spins in this picture: An individual bubble 
can be treated as a point-like object moving through a sea of spins. Since the exchange interaction, 
J, in the WC is so small, this is analogous the Nagaoka problem in that there is a strong tendency 
of to align the spins in a region (a ferromagnetic polaron) around each bubble so as to minimize 
the bubble zero-point kinetic energy. At T = 0, the size of this polar on [t^ is 

Lp ~ [nMeJ/h'^]-^/^ (33) 

where Mb is the bare effective mass of the bubble (given by an expression similar to that in Eq. 
I32() . Eq. ESI has been derived under an assumption that Lp ^ Lb, or that the exchange interaction 
J is small enough. In this limit, the bubble effective mass is strongly renormalized, most probably 
[7^ by an exponential factor Mp — > Mp exp l+AnLj:,] . 

However, if Lp is large, even at rather small bubble concentration the polaronic clouds begin 
to overlap, implying that there should be a transition to a globally ferromagnetic state. As in the 
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Nagaoka problem, the characteristic energy (temperature) scale of the ferromagnetism can be large 
compared to J. 

B. Delocalization of WC bubbles in the presence of slight disorder. 

If p <C h/e^, the weak localization corrections to the conductance are negligible on practical 
length and energy scales. However, WC bubbles are bigger and heavier than electrons. The 
question arises whether they are localized by weak forzen disorder. The purpose of this subsection 
is narrow: we will show that for small enough p, the WC bubbles are not localized by weak disorder. 
The reason for this is that the WC bubble can quantum melt and re-cryslallise. 

To a first approximation, the bubbles are pinned at favorable sites of the disorder potential, 
which we label with a discrete index, j. Since the Fermi liquid degrees of freedom are "fast," we 
can integrate them out to obtain an effective hamiltonian for the bubbles. In this case a relevant 
effective bubble Hamiltonian to consider has a form 

H-ff = J2 e^a^al + \^fala,a]a, + J] [ajaj + H.C.] + . . . (34) 

Here, creates a bubble at position R^. Since by assumption the disorder is weak and the bubbles 
are large, the optimal bubble at each site are all roughly the same size. Lb', the major effect of 
disorder is to pin their center of mass motion. We assume that the energies are randomly 
distributed in some interval of energies G [— VF, VF]. According to Eq. |S1 at long distances, 
V^-^j ~ |R*-Rj|-^ In the mean-field groundstate, the sites with < e are occupied by a single 
bubble (The repulsion between two bubbles on the same site is effectively infinite.), while the sites 
with Ci > e are "empty" , which means that the system at these sites is in a liquid state. 

The last term in Eq. which gives the bubbles quantumdynamics describes processes of 
transfer of bubbles from an "occupied" state at site "i" to an "empty" state at "j", or , in other 
words, processes where WC bubble melts at the site i and coherently recrystallizes at site "j". 
While the WC itself does not need to propagate long distances in this process, the density deficit 
associated with the bubble still needs to quantum diffuse through the FL. As a result, as shown in 



Appendix IbI using the methods of Ref. |4J], as |Rj — Rj| — > oo 



'''^ ~ ^^^^ 

where /3 ~ {AnLB)'^{pe^ /h) is a constant which vanishes as p ^ 0, and B is an exponentially small 
but non- vanishing coefficient. 

The crucial feature of this problem, which distinguishes it from the analogous problem of An- 
derson localization of electrons, is that the hopping matrix tij is long ranged. It is well known [3| 
that for (3 < D, Anderson localization does not occur. One can see this from the following: Let us 
consider an interval of energy je^ — e| < 5 <^W. A typical distance between sites with energies in 
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this interval is ^ {5u) where is a density of states per unit energy. Thus the ratio of the 
typical hopping matrix element to the energy difference between two such sites 

U,/5 ~ (36) 

diverges as 5 ^ 0, provided (i > D. Thus, if p is small enough, the bubbles are always be able to 
tunnel between nearly degenerate sites and hence are delocalized. 

The conclusion that the bubbles are delocalized for small enough p is apparently unavoidable. 
Taking the Hamiltonian in Eqs. OUandlSSla't face value, the present analysis suggests the existence 
of a bubble localization-delocalization transition at /> ~ /i/e^AnL^. However, more complicated 
processes, involving the coherent motion of many bubbles, could affect this conclusion. Finally we 
would like to mention that our calculations are controlled only in the limit when AnL^ ^ 1, where 
the amplitude B in Eq. ESI is exponentially small (See Eq. IB7|) . and consequently the mobility of 
the bubbles at T = is correspondingly small. 

VI. HYDRODYNAMICS OF MORE COMPLEX MICRO-EMULSION PHASES 

In the preceding sections, we have focussed primarily on electronic microemulsion liquid phases 
with no spontaneously broken symmetries. At low temperatures, the same physics leads to a 
variety of phases with different patterns of symmetry breaking. Because these phases are made 
out of clusters of electrons, rather than out of point-like electron quasiparticles, they are more 
analogous to the phases of complex classical fluids than to simple liquids. 

Microemulsion phases with different symmetries have different hydrodynamics. In this section, 
we will discuss a few simple aspects of the problem. Hydrodynamics deals with a set of "slow" 
variables, whose dynamics are slower than all microscopic relaxation rates. It can be used in 
an asymptotic sense since the variations of hydrodynamic variables are arbitrarily slow in the 
long wave-length limit. Examples of such hydrodynamic quantities are the temperature, average 
velocity, density, and the director in a nematic or smectic liquid crystal. The resulting sequence of 
hydrodynamic equations are the conservation laws plus the linear response equations connecting 
the fluxes with gradients of the parameters. These equations are asymptotically exact at times 
long compared to the characteristic relaxation time. However the term hydrodynamics can also 
be used in a slightly broader sense in systems with a substantial separation of time scales to 
include the dynamics of any slow variables on intermediate time scales. Below we will use the term 
hydrodynamics in both senses. We merely touch on the subject, to illustrate the richness of the 
possibilities. 
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1. "Supersolid" hydrodynamics 

Let us consider a case in which a dilute Hquid of FL droplets is dissolved in a WC host. This 
state has an unusual combination of liquid-like and solid-like properties. On the one hand, the 
system is characterized by a finite shear modulus. On the other hand, because of existence of the FL 
droplets, the state can flow, even in the presence of dilute pinning centers. To make the discussion 
as simple as possible, we will consider the case in which the electron spin degrees of freedom can 
be ignored; this pertains when the electrons are fully polarized by a magnetic field parallel to the 
film. The remaining slow variables are related to the density and velocity of thermally excited 
plasmons in the WC, the bubble velocity and density, and the elastic deformations of the WC, 
itself. The dynamics of this system is characterized by bubble-bubble Tfcf,, plasmon-plasmon Tpi^pi 
and bubble-plasmon Tb^pi scattering mean free times. Among them, a special role is played by 
umklapp-like processes in which a momentum equal to h times the WC reciprocal lattice vector is 
transfered from the of the bubble or phonon fluid to the WC lattice. 

At times longer than the umklap mean free time tjj, this problem is similar to any other crystal 
with a finite concentration of vacancies or interstitials: the hydrodynamic modes consist of under- 
damped longitudinal and transverse plasmon modes (where in the Coulomb and dipolar cases, the 
longitudinal modes have a Vk and k dispersions respectively) and a diffusive density mode. Here 
k is the wave vector of the plasmon. At sufficiently low temperatures, umklap processes become 
ineffective. In this case, on intermediate times tu ^ t ^ Tpi,phhTpi,bTTbb scattering processes con- 
serve the total quasi-momentum in the system of phonons+bubbles. In this case the hydrodynamic 
equations have the same form as in the case of conventional crystals where tjj ^ Tpi^pi 28]. Among 
other things, this can lead to interesting phenomena resembling second sound in the system of 
bubbles+phonons . 

At low enough temperatures, when the bubble scattering rate is smaller than their energy, 
bubble motion on some distances has a quantum coherent character. However, as long as the 
bubbles are not degenerate, the quantum nature of their motion manifests itself only in the values 
of the transport coefficients without changing of the form of the hydrodynamics. 

At even lower temperatures, new behavior emerges which depends on the bubble statistics. 
For instance, if the bubbles are bosonic, they bose-condense at low temperature. In this case the 
system of hydrodynamic equations is equivalent to the supersolid equations 123^. If the bubbles 
are fermionic, the system is analogous to a CDW in which a portion of the Fermi surface remains 
ungapped. One of the most important consequences of the existence of a bubble fluid is that, even 
at T = 0, the system can flow in the presence of obstacles and, therefore, it can conduct. This is 
very different from the case of a conventional WC, which is pinned by arbitrary small disorder. 
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2. Crystallization waves in the stripe phase. 



To illustrate a different aspect of the hydrodynamics of microemulsion phases, we consider a 
stripe phase. Since the symmetry of the stripe phase is the same as that of a smectic liquid crystal, 
at long enough times, the hydrodynamic equations should be identical to the equations describing 
the hydrodynamics of a smectic. 

On the other hand, on intermediate times, a reduced description of the system is also possible. 
One has to solve the Navier-Stokes equations for a FL and the elasticity theory for WC stripes 
and supplement these equations with boundary conditions at the WC-FL interface. The nature of 
these conditions at T = (whether the surface is quantum rough or quantum smooth) has been 
discussed in the context of the theory of '^He 01. In the quantum rough case, in the ground state 
there are steps on the boundary and it moves in response to the infinitely small jump in chemical 
potential between the FL and WC. In the quantum smooth case, there are no steps at T = 0, so 
the boundary cannot move until the difference in the chemical potentials reaches a threshold value. 
At finite temperature there are always steps on the surface. Therefore, it can move, and there is 
no sharp distinction between the quantum rough and quantum smooth cases. 

The long wave length dynamics of the system turn out to be independent of the character of 
the interface, including the concentration of the steps on the surface The reason is that the 
densities of the crystal and the liquid are different and as the surface moves, the additional mass is 
released into the liquid. Thus the effective mass associated with the long wave length motion of the 
surface is determined by the mass of the moving liquid rather than by the contribution associated 
with the motion of the surface itself. 

More generally, the fact that the two different components of the electronic microemulsion 
can inter-convert leads to qualitatively new modes, not present in classical liquid crystals. We 
illustrate this by considering, for the stripe phase, the simple situation when the crystallization 
wave propagates perpendicular to the stripes. Let us introduce the deviation of the position of 
the WC-FL boundary from its equilibrium position h{t). Then, according to Eq. I^Jthe potential 
energy per unit length of stripe is Vp^^ ~ ah'^ fwc / Lq. The kinetic energy is associated with the 
uniform flow of the liquid between two stripes. The flux density of mass released into the liquid 
by the motion of a boundary is J = mAnh{t) = mnU, where U is the velocity of the liquid. As a 
result, at fwc 1 the kinetic energy of FL is = mnLQf^(j{/\n)'^h?{t) and hence the frequency of 
oscillations is 




(37) 
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VII. DRAG EFFECT IN DOUBLE ELECTRON LAYERS 



In this section we discuss the more comphcated situation in which there are two 2DEG's sepa- 
rated by an effectively impenetrable barrier, but which interact with each other by the Coulomb 
interaction. The phase diagram of the two layer system is considerably richer than that of the sin- 
gle layer, even at mean-field level, due to the possibility of different but strongly coupled patterns 
of Coulomb frustrated phase separation in each of the layers. We will discuss a partial solution of 
this problem in a forthcoming paper j?^. 

In the present section, we will study the two layer system, because it is the geometry that 
permits a class of very revealing experiments of the "drag" resistance, which allow us to obtain 
rather direct information about the character of correlations at intermediate length scales. The 
drag resistance, pD = Vi/12-, is the ratio between the longitudinal voltage induced in the upper 
layer Vi by a current in the lower layer I2. To be concrete, in this article we will discuss only the 
case in which one of the layers is in a bubble liquid phase. Since the mobility of the WC bubbles in 
the lower layer is finite (even at T = 0) , they move with a velocity proportional to Ii . Images of the 
bubbles scatter electrons in the upper layer, thus transfereing momentum to the electron system 
of that layer. Generally, at small temperatures this contribution to the drag resistance is much 
bigger than that in the Fermi liquid. Its value depends on the ratio d/Ls- Specifically, if d is the 
spacing between layers, the drag effect is large when the characteristic size, of the mobile charged 
objects in the lower layer is greater than or equal to d, while relatively little contribution to the 
drag resistance comes from charge motions which are correlated on length scales small compared 
to d. Thus, the magnitude of pD is an indirect measure of the characteristic size and magnitude 
of the charge fluctuations in the system. 



A. Drag between two Fermi liquids. 



As a warmup, we briefly review some of the results concerning the drag between two pure Fermi 
liquid electron layers. In this case the electrons can be described by the Fermi liquid equations 

9ek54^'^^(r) 5e9n^;^ 1,2. ^ .^^s 
+ 7^ = ^st (^k) (38) 



^k 

,(1,2) 



9k dr Ok dr 



ek + E^kk'4' (39) 



p 



where ' are quasi-particles distribution functions in the upper and lower layer respectively, 

(12) 12 

' are electron energies in the layers, ■^k.k' ^-I'e Landau functions, and 1^1 are scattering integrals 
describing inter and intra layer electron scattering. 

Let us assume that the electron liquid moves in the lower layer with a velocity U. Then the 
drag resistance is the ratio between the rate of the change of the momentum in the upper layer 
and the current j = eNU in the flrst layer. 
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According to Eg 1391 at finite value of U tlie electron spectrum of quasi-particles in the upper 
layer ek is anisotropic. Nevetherless at zero temperature the drag resistance is zero. Indeed, if the 
quasi-particle distribution function is an arbitrary function of the total quasi-particle energy ek 

J = E = (40) 

k 

The drag resistance at finite temperature between two Fermi liquid layers has been analyzed in 
detail in Ref. [4^. So long as the elastic mean-free path, i ^ d, the drag is given by its value in 
the absence of disorder , 

PD ~ {kFd)-\kF/qTF)\T/EFf, (41) 

where qtf = 2^^/'^rskp is the inverse of the Thomas Fermi screening length. Thus, pD increases 
quadratically with temperature but is small both because the density of current carrying excitations 
is small (~ T/Ep) and because the length scale which characterizes the associated fluctuations in 
the charge density is small compared to the spacing between layers (assuming kpd S> 1). Density 
fluctuations in the Fermi liquid are enhanced by disorder, with the result that (in addition to a 
term which is logarithmic in T), pp), although still small, is enhanced by a factor of {d/i)"^ in the 
limit i<^d. 

The fact that p/j^OasT^Oisa property of the Fermi liquid and is not a generic property 
of all liquids. To see this, consider the drag between two superfluids. In this case, the expression 
for the super-current of the second layer 

Js2 = A^s2U,2 + A^sl2U,i (42) 

has a part which is proportional the the superfluid velocity in the first layer Usi Here \]sa is 
the superfluid velocity in the a*^ layer. The cross-coefficient A'^si2 in Eq- EH has a nonzero value at 
r = 0. Thus the system of two super-fluids exhibits non-vanishing drag at T = 0. 



B. Drag between two liquids, one of which is in the bubble phase. 

Let us consider the case when the electron densities in the layers are different. Specifically, let 
us consider the case in which the "upper" layer has a large enough electron density (small enough 
Tg) that it forms a uniform FL, and thus plays the same role as the metallic gate in a MOSFET. 
To begin with, let us consider the case in which the "lower" layer forms a bubble liquid in which 
the minority phase (bubbles) is a WC. As shown in Sec. IIVI in this case there are two types 
of the current carriers: electrons and the crystalline bubbles. At zero temperature the electron 
Fermi-liquid contribution to the drag resistance is always small and vanishes as T — > 0. 

Let us consider the bubble phase in the temperature regime in which of the upper layer is 
greater than the spacing between bubbles, riB^le ^ 1. In this case, transport in the upper layer 
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can be treated using the Boltzman equation. Each Wigner crystal bubble in the lower layer has 
an image of opposite charge in the upper layer. Electrons in the upper layer scatter on the image 
potential of the bubble, so the relative motion of the bubbles in the lower layer leads to a non-zero 
momentum transfer between the layers, even at T = 0. Specifically, in the absence of disorder 
pinning, the bubbles drift with the same mean velocity as the rest of the fluid in the lower layer. 
Thus, 

PD = — o"B (43) 

n 

where as ~ max[d, pb] is the cross-section for scattering the electrons in the upper layer on the 
image potential from the bubbles in the lower layer. In the presence of disorder, the large, heavy 
bubbles will be less mobile than the rest of the fluid in the lower layer, so po in Eq. 231 will be 
reduced by a factor of the relative mobility of the bubbles. However, we have shown that bubble 
localization does not occur, for weak disorder so, although pjy may be very small, even as T ^ 
it remains finite. 

Other cases, such as the hydrodynamic case in which l^e ^ Lb, have been discussed in Ref. [s^- 
The qualitatively important point is that in all cases, pD oc ns, so pD is directly sensitive to the 
characteristic T and ffn-dependences of nB{T,H) which originate from the Pomeranchuk effect. 



The d-dependence of the effect has also been discussed previously by us in Ref. l22l . At sfficiently 
small temperatures T the value of /jd in microemulsions is always larger than its Fermi liquid value 
and increases as the parameter Lb/cI increases. 

VIII. A COMPARISON WITH EXPERIMENTAL RESULTS 

In this section, we return to the list of experimental signatures of non-FL behavior that we 
discussed in Sec. ^ and compare them to the behaviors expected on the basis of the theory of 
electronic microemulsion phases. (Recall, the electron liquid Uj , Id, p, \A I sL Ifi , pi ™ Si 

MOSFET's spans a range of = 5 - 20 and the hole hquid[ll|, UJ, UJ, UJ, uA uA UJ in p-GaAs 
heteroj unction spans = 10 — 40.) We treat the experiments in the same sequence they were 
presented in Sec. ^ In the Section IVIII CI we discuss problems we see with applying some of the 
alternative, FL based theories to the interpretation of these same experiments. 

We have established that, in the absence of disorder and for a range of in the neighborhood 
of a critical value, electronic microemulsion phases must occur. The most unambiguous exper- 
imental evidence of such phases would come from the detection of an accompanying pattern of 
spontaneous spatial symmetry breaking. There has been, to date, no direct and unambiguous 
experimental confirmation of these phases at H± = 0, although related effects have been observed 
in the ultraclean 2DEGs in the presence of a large H± js^. However, quenched disorder is always 
a "relevant perturbation," which makes macroscopic observation of symmetry breaking difficult; 




FIG. 7: Schematic representation of the expected temperature dependence of the resistivity of an electronic 
microemulsion phase in which the FL component is the majority phase. The upper curve corresponds to 
H\\ > H* while the lower curve corresponds to — 0. See discussion in Sec. Ill C II 



probably, spontaneous symmetry breaking has not been detected in the current generation of large 
rg devices because their mobility is still insufficiently high. 

It is also worth acknowledging that the microscopic aspects of the theory are not, at present, 
under good quantitative control. Not only does the characteristic length scale, Lq, depend sensi- 
tively (exponentially) on unknown microscopic parameters such as the surface tension, but even 
the range of over which the microemusion phases occur is uncertain. The numerical calculations 
that lead to the generally quoted value of r* = 38 [s^ do not take into account the fact of the 
existence of the microemulsion phases and the enhanced stability of the WC at finite T due to the 
spin-entropy {i.e. the Pomeranchuk effect). Therefore, even this basic issue must, at present, be 
treated phenomenologically. 

Below we discuss a set of experimental data which indirectly indicate the existence of FL - 
WC phase coexistence in the strongly correlated 2DEG and 2DHG with large r^. The theoretical 
interpretation of these data are based mainly on the Pomeranchuk effect. As discussed in Section 
III C 11 this leads very generally to the qualitative behavior of the resistance as a function of T and 
i^ll shown schematically in Figs. [7|and|Hl 
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FIG. 8: Schematic representation of the H\\- dependence of the resistivity of strongly correlated electron gas 
under the same circumstances as in Fig. |7| 



A. Experiments on single layers 



a. The existence of an apparent metal-insulator phase transition: For zero disorder, there are 
several distinct phase transitions between different electronic microemulsion phases which occur as 
a function of r^. However, even weak disorder rounds the symmetry breaking transitions. However, 
a sharp "percolation" type transition between conducting and insulating states still occurs when 
the WC droplets overlap and block the electron transport through jthe FL. This, in principle, can 
qualitatively explain the metal-insulator transition observed in 



ough the FL. This, in princi ple, 



The Pomeranchuk effect implies that the critical electron concentration for the metal-insulator 
transition is a decreasing function of i^y, as is observed in experiment. 

b. T dependence of p{T,H\^\ = 0) at large rs and p <C h/e^: The significant increase of 
the resistance of metallic samples ( p(T) < C h/e'^ 



increasing temperature 



ic samples {p(T)< 

mmm 
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at relatively low temperatures (T < Ep) with 
can be interpreted as being a consequence 



of the Pomeranchuk effect: As the fraction of WC increases linearly with increasing T, so does p. 
As discussed in Sec. Ill El at T ~ Ep, the WC melts and a crossover occurs to a semiquantum 
regime in which the viscosity Eq. 1221 and consequently, the resistance fall in proportion to l/T 
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or faster (Eq. I23|) . Thus, theoretical expectations are consistent with the fact that p{T) reaches a 
maximum and then decreases when T > Ep. 

c. dependence of p{T ~ 0,-ff||) at large and p <^ h/e^: The Pomeranchuk effect impHes 
that the low temperature resistance of the 2DEG increases linearly with H\\ for i^y < Ep and sat- 
urates when > Ep^ where complete spin polarization is achieved. All this is in good qualitative 
agreement with the experimental observations in Refs. llSH . For small i^ii 
there should be a scaling relation between the T and //|| dependences implied by Eq. 121 which has 
not, to the best of our knowledge, been critically tested in experiment. 

d. The quenching of the T dependence of p{T,H\\) for > H* at large and p <^ h/e^: Close 
enough to the FL side of the microemulsion phases, the FL portion of the system can percolate even 
when the fraction of WC has been enhanced by the application of > . Thus, far enough from 
the "metal-insulator transition," we expect a range of densities in which the resistance increases 
strongly as a function of //|| , but still has a metallic value at high field: fi/e^ ^ p{H*\) ^ ~ '^)- 
However, once all the spins are polarized (for H\\ > -ff|p, the Pomeranchuk effect is quenched, and 
consequently so is its contribution to the T-dependence of p{T). Thus, in our opinion, the most 
direct confirmation of the central role played by the Pomeranchuk effect comes from a comparison 
of the temperature slopes of the resistance with and without the magnetic field. The T dependence 
of p is found to be 1-2 orders of magnitude weaker for ffy > H* than for ffy = and can even 
change sign. 0, Q| Such a dramatic effect strongly suggests that H\\ completely quenches the 
dominant contribution to the T dependence. 

e. '^Metallic" T dependence of p{T,H\^\ = 0) at large and p > h/e^: The conventional 
picture of the disorder driven metal-insulator transition implies that samples with large resistivity, 
p » /i/e^, at low temperatures should exhibit hopping conductivity, so p is expected to increase 
strongly as the temperature decreases. However, there are cases (See for example Fig.l in Ref. 
[l^ l in which highly resistive 2DEGs (p ~ (2 — 3)/i/e^) with large violate this expectation, in 
that p is an increasing function of increasing T. If what is changing with T is not a scattering 
rate, but the fraction of the system that is highly conducting, then there is no maximum value of 



the resistivity for which such metal 
wider implications for the so-called 
materials, for which related 



9C, 



ic T dependence is possible. This observation may have much 



95l | "bad metal" behavior observed in many highly correlated 



9ll ] intuitive explanations have been proposed. 



f. Oscillations of p{H±) at T > Ep ■ The existence of the oscillations of the resistance p{H± 



as a function of H±^ at ilp S> T ^ Ep is an experimental fact which drastically contradicts 
FL theory. At these temperatures, the WC has melted and the system is in the semiquantum 
regime. Although we do not have a theory explaining this phenomenon, we suspect that it reflects 
directly the strongly correlated nature of the electronic system: Since the kinetic energy in the 
electronic system is still much smaller than the potential, on short spacial scales the liquid may 
have correlations which are similar to those of a crystal. Thus, in the interval of magnetic flelds 
where the magnetic length is comparable to the interelectron distance, short range quaisiper iodic 
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order may account for a few oscillations of p{H±). 



B. Experiments on double-layers 

We now turn to the anomalies in the measured drag resistance po in p-GaAs double layers 
with small resistances p ~ (1/10 — l/20)h/e'^ and large r^. 

a. The magnitude of p£): The fact that pz) is 1-2 order of magnitude larger than expected on 
the basis of Fermi liquid theory is in qualitative agreement with our theory. Specifically, where the 
FL is the majority phase in the lower layer, the fact that the expected size of the WC domains ( 
Eq. 0} is Lq > d, where d is the spacing between layers, means that none of the familiar factors 
responsible for the small size oi pD in a FL (Eq. I41() apply. Put another way, the large magnitude 
oi pd can be viewed as direct confirmation of the existence of dynamical charge inhomogeneities 
of substantial magnitude and on length scales greater than or equal to d. 

b, c, and d. T and H\\ dependence of pd{T,H\\), and its similarity to p[T,H\\): In addition to 
its large magnitude, the most striking feature of the measured /9_d in large double layer devices is 
that the T and especially the H\\ dependencies looks qualitatively similar to the resistances of the 
individual layers p{T,H\\). Most notable is the fact that, again, a strong magnetic field, H\\ > H* 
substantially suppresses the temperature dependence of Pd(T) in qualitative agreement with 
the assumption that the dominant T dependence is a consequence of the Pomeranchuk effect. 
More generally, the qualitatively similarity between p and pn speaks in favor of the assumption 
that properties of both p{T,H\i) and p£){T,Hii) are connected with phase separation. 



C. Alternative theories of experiments in large rg devices 
Attempts have also been made (Refs. 0,13,1411) 

to explain the observed behaviors by extrap- 



olating known perturbative expressions which are valid in the limit of small r<j and p <^ h/e^, to 
the regime of large r^. In this section, we point out some of the difficulties with these FL based 
explanations. Before getting to specifics, however, it is worth noting from the beginning that there 
is a general problem with applying these theories to the experiments discussed above: Perturba- 
tive corrections are, by their nature, small compared with the resistance itself, while the observed 
variations of p with and T are large - up to a factor of 10. Moreover, in all cases, to explain the 
experimental data, the FL results are being extrapolated far beyond their natural range of validity 
— specifically, is much greater than 1. 

In the weakly interacting limit, ^ 1, and at kp£ ^ 1, calculations of the resistance of the 
2DEG can be carried out in a controlled fashion. To lowest order in powers of (kpi)^^ ^ 1, the 
transport properties of the system can be described by the Boltzmann kinetic equation and the 
Drude formula for the resistance, which at low T is determined by the impurity scattering. The T 
and H\\ dependences of the resistance p(T, H\\) are determined by the corresponding dependencies 
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of the density of states at the Fermi level, the Fermi velocity, and the scattering cross-section, and 
typically p{T,H^^) decreases weakly with T and H^^ because the electron-impurity cross-section is 
a decreasing function of the electron energy. Even the sign of the effect is opposite to what is 
observed in large devices. 

If p <^ h/e^, and at moderately low temperatures, both weak localization 0, and the 
electron interaction interference corrections [49:] to the classical Drude resistance are parametri- 
cally small. The physics of the interactional corrections is the following: Impurities in a metal 
create Friedel oscillations of the electron density. Due to the electron-electron interactions, the 
quasiparticles in the metal are scattered not only from the impurity but also from the modulations 
of the electron density. The interference between multiple electron scattering on impurities and 
on the Friedel oscillations of the electron density give rise to the non-analytical in T and H\\ cor- 
rections to the Drude resistance. At high enough temperatures when Lt = vp/T ^ / the former 
effect is reduced to an interference between two scattering amplitudes: The electron scatters from 
a single short range impurity from the Friedel oscillations created by the same impurity. The 
interference corrections to the Drude formula in this so called "ballistic regime" were considered in 



37 



3S 



3£, 



40|. The result is that in an intermediate interval of temperatures Ep ^ T ^ the 



T- dependent correction to the resistivity {p{T) — p(0)) ~ CT is linear in T and C ~ Fo(l + 15i*b) 
[4^. Here Fq is a Fermi liquid parameter which, for r<j <C 1, is negative with a magnitude pro- 
portional to Vg. Thus, at small r^, when Fq is small, C < and p{T) decreases as T increases. 
However, at large enough r^, when Fq < —1/15, the coefficient C changes sign, and p{T) becomes 
an increasing function of T, which is in agreement with the experiments. In principle, at ~ 1 and 
-ff||,r ~ Ep the corrections could be of order one. It is hard, however, to imagine that the Friedel 
oscillations of the electron density induced by an impurity can make a much larger contribution to 
the scattering cross-section than the scattering from the impurity itself. Therefore, we think that 
this process (even when extrapolated to large r^) cannot explain changes in the resistivity by more 
than a factor of two, as are observed experimentally. 

Another aspect of the problem is that the calculations of Ref. j^^l are valid only for T and 
<C Ep, where one can neglect the energy dependence of the electron scattering cross-section. 
In addition, for r<j <C 1, where the theory is well controlled, at i/y > Ep, or at T ~ Ep, the main 
contribution to the magnetoresistance p{H) — p{0) comes not from Friedel oscillations considered 
in 0Lbut rather from the ffy dependence of the Drude part of the resistance. At last the theory 

3, S S, 

ot explain why the T-dependence of the resistance p(T) is completely 



Ref. 



quenched at Hu > H*. Therefore we think that, though the results of the perturbation theory 
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3c 



40| may be relevant to experiments performed on samples with relatively small and 



at relatively small values of i/y and T, they can not explain results of experiments at ^ 1 where 
almost all effects are larger than order one. 

An attempt at explaining result, on large samples was undertaken in Ref. using the 

classical formulas for the resistivity and extending them to the regime ^ 1 and kpi ^ 1. There 
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are several aspects of this extrapolation that are troubling. Firstly, as a function of r^, both dp/dT 
and dp/dH\\ change sign, so the extrapolation cannot be said to be featureless. Moreover, the 
origin of the sign change can be traced to the fact that at large r^, the Thomas Fermi screening 
length is parametrically smaller than the spacing between electrons, kxF ~ fskp ^ kp- Screening 
lengths less than the distance between electrons are clearly unphysical. Moreover, the calculation 
treats the impurity scattering in Born approximation, which is also unjustified in this case. On 
a phenomenological level, the extreme degree to which the magnitude of dp/dT is suppressed by 
large H\\ in both p-GaAs and Si MOSFETs requires an unreasonable amount of fine tuning of 
parameters. 



IX. CONCLUSION 



In this section we discuss several important unsolved problems related to Coulomb frustrated 
phase separation. 



A. Unsolved problems concerning the equilibrium phase diagram 

a As discussed in Sec. quantum fluctuations necessarily melt a bubble crystal when the 
spacing between bubbles is sufficiently large. Thus, the mean-field Lifshitz transition between the 
uniform fluid and the bubble crystal is destroyed by quantum fluctuations; the bubble crystal 
phase is only stable when the lattice constant is smaller than a finite critical value. However, any 
transition from a liquid phase to a crystal with finite period is expected to be first order, in conflict 
with the theorem forbidding the existence of first order transitions. This is a conceptual problem 
that we have failed to resolve. 

A related unresolved problem is the quantum (T = 0) properties of the melted WC bubble 
phase. 

b. The properties of the interface between the FL and the WC are largely unknown. The 
question of whether the interface is quantum smooth or quantum rough is open. In addition, the 
characteristic periods of the two phases are different - 271 /Ikp for the FL and the lattice spacing 
for the WC, opening the possibility of interesting effects of incommensurability. The interfaces 
may even have phase transitions of their own. 



B. Significance of quenched disorder 



The strongly correlated liquid in the presence of an external scattering potential is a largely 
unsolved problem. There is a theorem forbidding first order phase transitions in 2D with arbitrarily 
small disorder Sj]- In the case of weak disorder, though, the characteristic length at which the first 



order phase transition is destroyed can be exponentially large. According to |52l |. whether the first 
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order phase transition is converted to a second order one, or to a cross-over depends on whether 
the two phases have different symmetries. However, long range crystahine order is impossible in 
the presence of quenched disorder in dimensionalities less than four [s^; the result is probably, a 
glassy state. 

One can see the glassy aspects of the problem even on the mean field level if one assumes that 
the coefficients in Eq. |21 are smooth functions of the coordinates. So long as Lq ^ ^, a coarse- 
grained description of the thermodynamic state of the system can be obtained by associating an 
"up" psuedo-spin with each local position which is in the WC state, and a "down" with each FL 
position. The problem of finding the mean-field ground-state then corresponds to the same problem 
in an effective random-field Ising model, albeit one with power-law interactions. While there is 
little that is actually known about the random field Ising model with longer range interactions, it 
is certainly likely that, in common with the short-range model, it forms a glass. 

Implications for localization theory in 2D conductors: An important open issue concerns the 
implications of Coulomb frustrated phase separation for the theory of localization in the presence 
of weak disorder. It is believed that at T = 0, all electron wave functions in disordered systems 
are localized. This statement is based on the logarithmic divergence of the interference corrections 
to the Dn,de conducUvityHQ Q and its interpretation via the renonnalization gtonp In 
turn, these results are intimately connected to the fact that, in the presence of an elastic scattering 
potential, the number of quasiparticles with a given energy is conserved. 

In the case when there is a finite concentration of WC droplets in the FL, the status of this state- 
ment is unclear. As discussed in Sec. due to quantum fluctuations, the bubbles can quantum 
melt and recrystalize. Thus the number of electrons in the Fermi liquid is not conserved because 
they can join the WC fraction of the system. Consequences of this fact for the 2D localization 
theory remain to be investigated. 



C. Mesocopic properties 



An important area of research involves transport through metallic "point contacts." Up until 
now, however, they have been studied only in three regimes: the regime where electron interactions 



are not important |83], the Coulomb blockade regime, and the Kondo regime jS^] in which there 
is a single localized electronic state through which tunneling between the two metallic reservoirs 
occurs. However, there is another interesting regime when there is a relatively large "depletion 
region" between two metallic reservoirs with high electron densities. In the depletion region, there 
is a strongly correlated electron liquid with low electron density whose density can be varied by 
changing a gate voltage. In this case as a function of gate voltage, aspects of the electronic micro- 
emulsion phases can be directly probed on a mesoscopic scale. 
In this context we would like to mention several problems: 

a) The Pomeranchuk effect can produce significant T and Hh dependencies to the resistivity 




FIG. 9: A schematic picture of a smooth gutter hke potential in a spht-gate device. The dots indicate 
electrons in the WC. The gray area corresponds to a melted FL bubble phase. 



through the depletion region. Suppose that at T = the system is in the FL state, but with 
the gate tuned so that the WC is close in energy. Then, as the temperature or the magnetic 
field increase, the resistivity can exhibit a sharp crossover from very conducting to very insulating 
values. This can mimic some of the behavior traditionally associated with the Kondo effect. 

b) An example of a mesoscopic geometry which is well suited to exposing electronic micro- 
emulsions is one in which the depletion region is produced by a smooth external potential which 
depends only on one coordinate. This "gutter-like" potential is shown qualitatively in Fig. |3 
Specifically, it is interesting to consider the nature of the transition from the insulating state when 
the entire gutter is filled with low density WC, to a conducting state in which the highest electron 
density portion of the system, which lies along the line of minimum potential at the center of the 
gutter, is a liquid. At T = 0, in the limit that the gutter is long, this is a one dimensional quantum 
phase transition. 

At mean-field level, the first transition that occurs with increasing electron density is to a phase 
with far separated FL bubbles; because the bubbles have relatively higher electron density than 
the surrounding WC, they are bound to the center of the gutter forming a ID bubble crystal. 
However, quantum fluctuations of the bubble's positions will always cause the bubble crystal to 
melt. The resulting bubble liquid phase is thus a conducting (fluid) phase. 

The bubble motion also produces changes in the spin-state of the surrounding WC. As discussed 
in Sec. O the FL bubbles are point like quantum-mechanical objects whose zero-point kinetic 
energy is lowest when the surrounding spin's are ferromagnetically polarized. Thus, as discussed 
in the neighborhood of Eq. 47, if the ratio of the bubble kinetic energy to the exchange energy 
in the WC is large enough, zero temperature kinetically stabilized ferromagnetism can onset with 
the onset of the bubble liquid phase. In this context, we note that there are a set of experiments 
Ref. 
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82 



83 



potential we have 
observed. In Ref. 



8J] in various mesoscopic devices, with a geometry similar to that of the gutter-like 
ust discussed, in which a set of anomalies referred to as the "0.7 effect" have been 
8l| , this anomaly was interpreted as evidence of the existence of ferromagnetism 
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FIG. 10: Schematic representation of the periodic osciUations of a WC droplet bound to a potential 
"mountain-top" which results in quantum pumping action, detectable in drag. 



in Id channels at large r^. 

c) Another interesting but still incompletely understood issue involves a particular type of 
quantum oscillations of the WC droplets of the sort illustrated in Fig. 1101 Here the black dots 
represent electrons in a crystal droplet which is bound at a preferred position by the disorder 
potential. Consider a process in which the bubble sloshes a bit to the right, then an appropriate 
sized sliver of the displaced crystal quantum melts, and an equivalent sliver of liquid coherently 
recrystallizes on the left. The result is that the bubble has returned to its original position. The 
value of the action corresponding to such a process is finite. An important feature of the process, 
however, is that there is net charge motion that accompanies it. 

In equilibrium, oscillations of the droplet are equally probable both to the right and to the left. 
However, the problem can be biased by application of a voltage. Then, this effect can be detected 
in drag: as the droplet in the lower plane executes this quantum motion, the image charge in the 
upper layer also acts as a pump of electrons from the left to the right. This gives additional support 
to the idea of a non-zero drag resistance in the limit T ^ 0. 



QC , in which the related 



Note added: After we completed this article, we became aware of a paper, 
problem of Coulomb frustrated phase separation in a bosonic system was studied by quantum 
Monte-Carlo methods. Here stripe and bubble crystal phases are clearly seen, although various 
subtleties, such as the hierarchical microemulsions at that we would expect between the stripe and 
bubble phases, have not yet been explored. 
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X.P.A. Gao, and D. Slialiar for useful discussions. 
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APPENDIX A: COULOMB FRUSTRATED PHASE SEPARATION 



In this appendix we present a more extended analysis of the results obtained in Refs. jlfll . 12(1 . l2ll | , 
give qualitative explanation of existence of the phase separation on the mesoscopic spacial scale 
(microemulsions), and derive Eqs. ^andl^l 



a. Qualitative Features 

A putative first order transition from the WC to the FL would occur when the free energy 
densities of the two uniform phases, Fwcin-) and FpLin-), cross, as in Fig. [l] In the case in which 
V{r) is long-range, i.e. it falls as or slower, F{n) must include the interaction with a uniform 
positively charged background which is assumed to be static and uniform. As a first step, we show 
that for n near Uc, neither a uniform phase, nor a macroscopically phase-separated combination 
of them, minimizes the free energy - rather some form of mesocscale phase separation always has 
lower free energy. 

To see the point most simply, consider the system with average density n = ric- We divide the 
system up into alternating stripes of WC and FL of width L ^ ^, where ^ is the width of the 
interface between a region of WC and FL. ( We consider stripes to be explicit, but other patterns, 
such as an array of bubbles of radius L, have similarly low free energies.) We also rearrange the 
density, so the WC stripe has density {ric — An) and the FL stripe has density (ric + An). The free 
energy density of this state can be expressed as a sum of different contributions 

F = Fun + Fps + F^ + Fir. (Al) 

The free energy density of the uniform density phase is F^n = Fwcin-c) = FpLinc)- The next term 
is what would have been the gain from phase separation, were there no long-range interactions 

Fps= (l/2)[F^i(nc + An)-FH/cK-An)] (A2) 
= -CiAn-FC2(An)2-h... 

Here Ci and ^2 are the first terms in an expansion in which . . . signifies higher order terms in powers 
of An. The third term in Eq. 1 is the short-distance contribution to the surface energy. 

Fa = a/L (A3) 

where a is the surface tension. The final term is the the long-range interaction in this particular 
trial state: 

Fir = (1/20) j dr dr'V{r - r')5n{r)5n{r'). (A4) 

where the difference between the local density and the average density 6n{r) = itAn depending 
on whether one is in a WC or FL portion of the system, is the total area of the system, and 
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/ dr (5n(r) = 0. Though the = 1 (Coulmb case) and x = 3 (the dipolar case) are the 

most interesting, physically, we start with the case 1 < x < 3 because the result is more easily 
established than in the limiting cases. 

The integral in Eq. IA4I has a different character depending on whether x is greater or less than 
x = 2. 

a) The case 1 < x < 2. In this case the integral is evaluated as: 

Fir = AQiAnfL'^-'' for 1< x < 2, (A5) 

where ^ is a positive number of order 1. Notice that the energy density diverges as L — > oo, 
reflecting the fact that macroscopic phase separation is impossible with long-range interactions. 
Consequently, if An is small (so that the higher order terms in Fps can be ignored) and L large 
compared to ^, 

F = Fun- CiAn + ao/L + AQ{AnfL^-''. (A6) 

This expression can now be minimized with respect to An with the result 

An = (Ci/2Ag)L-(2-); F = F^n- L'^ [{Cl/4AQ)L^-^ - a] . (A7) 

Notice that the long-range interaction produces something resembling a scale-dependent surface 
tension, which is negative for large enough L. Minimizing this expression with respect to L, we 
find the characteristic scale for nano phase-separation, 

L = Lo = [AAQa/{2 - x)Cf]^/^''-'l (A8) 

This expression is valid if a, and hence Lq is sufficiently large, so that An is small; if a is not large, 
the result depends on the higher order terms in Fpg, although this typically will not change the 
result qualitatively. 

b) The case, 2 < x < 3.' In this case both the long and short-distance contributions to the 
integral are singular, with the consequence that 

Fir = Q{Anf [aT^"^"^^ - AL-^''-'^^ + ...] for X > 2 (A9) 

The first term can be thought of as a positive contribution to Fpg, C2 — ^ C2 + AQ^~^^~'^\ while 
the latter can be thought of as a negative, but again scale dependent contribution to the surface 
tension, a ^ a — AQ{An)'^L^'^ . 

In this case, macroscopic phase separation is possible, and the density difference, An* between 
the two coexisting phases is determined by minimizing Fps{An) = Fps{An) + ^'Q(An)^^^~^ with 
respect to An - i.e. by the Maxwell construction. The free energy of the macroscopically phase 
separated state that results is then F„„ + Fps(An*) < However, taking L to be large but 

non- infinite, results in a state with free energy density: 

F = Fur^ + Fp,{An*) - L-^ \AQ{An''fL^-'' - a] (AlO) 
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which is manifestly lower than the free energy of the macroscopically phase separated phase for 
any L > [a / AQ{l^n^)'^Y^^'^~^\ The optimal domain size is 

L = Lo = [a/{2 - x)Ag(An*)2]i/(3-x)_ (All) 



The cases of x = 1 and x = 3 are discussed in detail in [l^, l21J, |2]| . They require a more subtle 
analysis in the sense that it is necessary to consider a form of phase separation in which the density 
varies within each domain. For instance, for x = 1, in the low density (WC) domain, the density 
varies as ric — An/r as a function of the distance, r, from the interface with the high density phase. 
When this is done, the result is to replace L^~^ in Eq. IA71 L^~^ in Eq. lAlOl bv log[L/^] and the 
expressions in Eq. IA8I bv 

L = Lo = ^ exp[iAQa/Cf + 1]. (A12) 

Here the short distance cut off, ^ is of order d in the dipolar case and of order ineterelectron 
distance in the Coulomb case. (This result is equivalent to Eq. 4.) 

It is worth noting that, in the Coulomb case, the micro-emulsion phase has lower free energy 
than the uniform phase, while for the dipolar case, the competition is between the macroscopically 
phase separated phase and the microemulsion. However, this is the only difference - the same 
logarithms enter the free energy expression for the microemulsion phase in both cases, and hence 
the character of the microemulsion is, surprisingly, similar for Coulomb and dipolar interactions. 

b. Derivation of Eqs. Q| and\^ 

So as to be able to obtain explicit results, we consider an effective free energy density describing 
a first order phase transition between two phases, a high density phase, in which a suitably defined 
order parameter (/)(r) = 1, and a low density phase in which <j){r) = —1: 

F = T[(j)] - C(n(r) - nc)</'(r) + -(n(r) - n) f dr'V{r - r')(n(r') - n). (A13) 



2 

Here T[(j)] is a local functional of (p with two degenerate minima corresponding to 0(r) = ±1. It also 
has a saddle point corresponding to a domain wall between the two phases of width ^ and energy 
per unit length o" > 0. The final terms determine the long wave-length profile of the density, n(r), 
where n is a fixed uniform background density, and Uc is the critical density. (Where the optimal 
density of the two phases is very different, higher than quadratic terms in powers of (n(r) — n) may 
be quantitatively but not qualitatively important.) 
Minimizing Eq. IA13I with respect to n(r) we get 



n 



(r)-n = C^"V(r), (A14) 



where y ^ is the operator inverse of V, which can be defined by Fourier transform form as (n(r) — 
^)k = C(l/^)0fc- Substituting this expression back into Eq. UT3l we get an expression for total 
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energy of a microemulsion phase 

Fme = Jdr |f [</>(r)] - C(n - n,)0(r) - (1/2) J dr' dr<P{r)V {r - r')ar'0(r')| • (A15) 

For 1 < X < 2, since — > cxd as |k| — > 0, the dominant effect of the long-range interactions 
is expressible as a non- local stiffness, so f = T and V{r) ~ C^Q^VkP"""'- For 2 < x < 3, the 
dominant contribution is local, resulting in a renormalization of F — > F = F — ((^^/2Vb)0^, and it 
is the subleading term which produces a non-local stiffness, V{r) ~ C^Q~^/kl^~^- The amusing 
result is that the form of the free energy depends only on |x — 2|, and in particular that V{r) ~ 1/r 
both for the Coulomb case, x = 1, and the dipolar case, X — o. 

In the limit when the widths of the interphase is smaller than the size of domains one can 
neglect the gradient of 0(r) inside the phases and we arrive to Eqs. ^andEl 

APPENDIX B: DERIVATION OF EQ.I^ 

Below we present a heuristic derivation of Eq. EHl but the same result can be derived more 
formally following the methods of Levitov and Shvtov(44l|. 

To understand the origin of the non-local tunneling represented by Eq. 1351 we begin by consid- 
ering a process in which a bubble of WC at site j quantum melts. (We will find that this process 
has a logarithmically divergent action - were the action finite, it would lead to a term proportional 
to [oj -|- Oj] in Eq. 1351 ) We imagine that there are two nearly degenerate states of the system, one 
in which there is a bubble at site j, and one with a slightly different energy, —ej, in which site j 
is covered with liquid. Since in equilibrium, the liquid has a slightly higher density than the WC, 
the transition between the two states must involve the motion of an amount of charge Q = eAnL^ 
in from "infinity," i.e. from a distance L^o which is the system size. The classical action for this 
process is thus the sum of three contributions: 

S = Smelt + Shydro + Sf,. (Bl) 

where, as in the usual theory of quantum nucleation. 

Smelt ~ crL%/vF (B2) 

is the action involved in melting the WC bubble locally, Shydro is the action associated with the 
motion of a charge Q in from infinity, and (with e = —ej) 

Se = ~ j dre (B3) 

records the cost of leaving the system in a higher energy metastable state. Since the charge 
redistribution takes place entirely in the FL and dominantly at long distances, Shydro depends only 
on the hydrodynamics of the FL. In Eq. IB31 the integral runs over the "imaginary" time interval 
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which characterizes the event. As we shah see, both Shydro and 5^ diverge as Loo — > oo. However, 
when apphed to the process of melting and recrystahization at a site a distance R away, the same 
analysis leads to a finite result, with Lqo replaced by R. 

A lower bound estimate of Shydro can be made by integrating over time the instantaneous excess 
energy associated with the charge excess, 

POO 

Shydro ~ / drEir). (B4) 
Jo 

In the short-range case (including the dipolar case), 

E{t) = {1/2)k J (iV[5n(r,t)]2 ~ K[An]'^[LB/L{t)]'^ 

where k is the compressibility of the FL, L{t) = V Dt oc \/t j p is the distance to which the excess 
charge can diffuse in time r, and D is the diffusion constant of the FL. In the Coulomb case, 

E{t) = Q^Lir), (B5) 

where L{t) = {e^ /h)p~^t. Since, in both cases, E{t) — > 1/r at long times, the resulting action is 
logarithmically divergent with system size, 

Shydro ~ /3 log [Loo /Lb] 

/3 OC /j[AnL|]2. (B6) 

An estimate of S^ ~ — ctq is obtained by noting that the characteristic time involved in the charge 
motion is tq = L^/L* in the dipolar case, and tq = L^{h/e'^)p in the Coulomb case. 

The generalization of this analysis that leads to Eq. [33 is now obvious. Replacing Loo in Shydro 
with Rij = |Rj — Rjl leads to the stated result. Clearly, 

B oc exp[-25',„e/t] (B7) 

is exponentially small but finite, as promised. 

There is an additional, potentially stronger (exponential) Rij dependence of tij inherited from 
S^ with e = ei — ej, which we have ignored in this discussion. However, since for the relevant pairs of 
sites, \ei — ej\ ~ J is small, and Rij is only large in proportion to the product \ei — ej\To ~ 6^ 

in the dipolar case and |ej — ej|To ~ 5^^^ in the Coulomb case; in neither case does S^ make an 
important contribution to the action. 

Our calculations are controlled only in the limit when AnL^ 1, the amplitude in Eq. IB 71 
is exponentially small, and consequently the mobility of the bubbles at T = is correspondingly 
small. 
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